science  for  a  changing  world 


In  cooperation  with  the  Texas  Department  of  Transportation 

Literature  Review  for  Texas  Department  of 
Transportation  Research  Project  0-4695: 
Guidance  for  Design  in  Areas  of  Extreme 
Bed-Load  Mobility,  Edwards  Plateau,  Texas 


Open-File  Report  2005-1234 
(TxDOT  Research  Report  0-4695-1) 


U.S.  Department  of  the  Interior 
U.S.  Geological  Survey 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1 .  REPORT  DATE  2.  REPORT  TYPE 

2005  N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Literature  Review  for  Texas  Department  of  Transportation  Research 
Project  0-4695:  Guidance  for  Design  in  Areas  of  Extreme  Bed-Load 
Mobility,  Edwards  Plateau,  Texas 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Department  of  the  Interior  U.S.  Geological  Survey  1849  C  Street, 

NW  Washington,  DC  20240 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

18.  NUMBER  19a.  NAME  OF 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  |J|J 

unclassified  unclassified  unclassified 

50 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Technical  Report  Documentation  Page 

1 .  Report  No.  2.  Government  Accession  No. 

FHWA/TX-04/0-4695- 1 

4.  Title  and  Subtitle 

LITERATURE  REVIEW  FOR  TEXAS  DEPARTMENT  OF 
TRANSPORTATION  RESEARCH  PROJECT  0-4695:  GUIDANCE  FOR 
DESIGN  IN  AREAS  OF  EXTREME  BED-LOAD  MOBILITY,  EDWARDS 
PLATEAU,  TEXAS 

7.  Author(s) 

Franklin  T.  Heitmuller,  William  H.  Asquith,  Xing  Fang,  David  B.  Thompson, 

Keh-Han  Wang 

9.  Performing  Organization  Name  and  Address 
U.S.  Geological  Survey 
Texas  Water  Science  Center 
8027  Exchange  Drive 
Austin,  Texas  78754 

12.  Sponsoring  Agency  Name  and  Address 

Texas  Department  of  Transportation  U.S.  Geological  Survey 

Research  and  Technology  Implementation  Office  Texas  Water  Science  Cer 
P.O.  Box  5080  8027  Exchange  Drive 

Austin,  TX  7873 1  Austin,  TX  78754 

15.  Supplementary  Notes 

Project  conducted  in  cooperation  with  the  Texas  Department  of  Transportation  and  the  Federal  Highway  Administration. 

16.  Abstract 

A  review  of  the  literature  addressing  sediment  transport  in  gravel-bed  river  systems  and  structures  designed  to  control  bed-load 
mobility  is  provided  as  part  of  Texas  Department  of  Transportation  research  project  0-4695:  Guidance  for  Design  in  Areas  of 
Extreme  Bed-Load  Mobility.  The  study  area  comprises  the  western  half  of  the  Edwards  Plateau  in  central  Texas.  Three  primary 
foci  of  the  literature  review  are  journal  articles,  edited  volumes,  and  government  publications.  Major  themes  within  the  body  of 
literature  include  deterministic  sediment  transport  theory  and  equations,  development  of  methods  to  measure  and  analyze  fluvial 
sediment,  applications  and  development  of  theory  in  natural  channels  and  flume  experiments,  and  recommendations  for  river 
management  and  structural  design.  The  literature  review  provides  an  outline  and  foundation  for  the  research  project  to 
characterize  extreme  bed-load  mobility  in  rivers  and  streams  across  the  study  area.  The  literature  review  also  provides  a  basis 
upon  which  potential  modifications  to  low- water  stream-crossing  design  in  the  study  area  can  be  made. 


17.  Key  Words  18.  Distribution  Statement 

Gravel-bed  rivers.  Bed  load.  Low-water  crossing,  Edwards  Plateau.  Texas  No  restrictions. 


19.  Security  Classif.  (of  report) 

20.  Security  Classif.  (of  this  page) 

21.  No.  of  pages 

22.  Price 

Unclassified 

Unclassified 

40 

3.  Recipient’s  Catalog  No. 

5.  Report  Date 
June  2005 

6.  Performing  Organization  Code 

8.  Performing  Organization  Report  No. 
and  OFR  2005-1234 

10.  Work  Unit  No.  (TRAIS) 

1 1 .  Contract  or  Grant  No. 

Project  0^-695 

13.  Type  of  Report  and  Period  Covered 
Technical  Report  Sept.  2003  to  June  2004 

iter  - 

14.  Sponsoring  Agency  Code 


Form  DOT  F  1700.7  (8-72) 


Reproduction  of  completed  page  authorized 


Literature  Review  for  Texas  Department  of 
Transportation  Research  Project  0-4695: 
Guidance  for  Design  in  Areas  of  Extreme 
Bed-Load  Mobility,  Edwards  Plateau,  Texas 


By  Franklin  T.  Heitmuller,  William  H.  Asquith,  Xing  Fang, 
David  B.  Thompson,  and  Keh-Han  Wang 


In  cooperation  withtheTexas  Department ofTransportation 


Open-File  Report  2005-1 234 
(TxDOT  Research  Report  0-4695-1 ) 


U.S.  Department  of  the  Interior 
U.S.  Geological  Survey 


U.S.  Department  of  the  Interior 

Gale  A.  Norton,  Secretary 

U.S.  Geological  Survey 

Charles  G.  Groat,  Director 


U.S.  Geological  Survey,  Reston,  Virginia:  2005 

For  sale  by  U.S.  Geological  Survey,  Information  Services 
Box  25286,  Denver  Federal  Center 
Denver,  CO  80225 


For  more  information  about  the  USGS  and  its  products: 

Telephone:  1 -888-ASK-USGS 

World  Wide  Web:  http://www.usgs.gov/ 

Any  use  of  trade,  product,  or  firm  names  in  this  publication  is  for  descriptive  purposes  only  and  does  not  imply 
endorsement  by  the  U.S.  Government. 

Although  this  report  is  in  the  public  domain,  permission  must  be  secured  from  the  individual  copyright  owners  to 
reproduce  any  copyrighted  materials  contained  within  this  report. 

The  contents  of  this  report  do  not  necessarily  reflect  the  official  view  or  policies  of  the  Texas  Department  of 
Transportation  (TxDOT).  This  report  does  not  constitute  a  standard,  specification,  or  regulation. 


Suggested  citation: 

Heitmuller,  F.T.,  Asquith,  W.H.,  Fang,  Xing,  Thompson,  D.B.,  and  Wang,  Keh-FHan,  2005,  Literature  review  for  Texas 
Department  of  T ransportation  Research  Project  0-4695 — Guidance  for  design  in  areas  of  extreme  bed-load  mobility, 
Edwards  Plateau,  Texas:  U.S.  Geological  Survey  Open-File  Report  2005-1234,  40  p. 


Contents 


Abstract . 1 

Introduction . 1 

Detailed  Literature  Review  of  Selected  References  . 3 

Literature  Reviews  by  Lamar  University  Researchers . 3 

SedimentTransport  Mechanism  in  Gravel-Bed  Rivers . 3 

Hydraulic  Geometry  of  Gravel-Bed  Rivers . 7 

Sediment  Transport  Equations  for  Gravel-Bed  Rivers . 7 

Laboratory  Experiments  for  Gravel  Transport . 8 

Literature  Reviews  by  Texas  Tech  University  Researchers . 10 

Literature  Reviews  by  University  of  Houston  Researchers . 10 

Debris-Flow  Managementfor  Low-Water  Crossing  Protection . 10 

Research  on  Debris-Flow  Countermeasures . 15 

Conclusions . 18 

Literature  Reviews  by  U.S.  Geological  Survey  Researchers . 19 

Summary . 22 

Bibliography . 22 

Annotated  Bibliography  by  Lamar  University  Researchers . 22 

Annotated  Bibliography  by  Texas  Tech  University  Researchers . 24 

Annotated  Bibliography  by  University  of  Houston  Researchers . 24 

Stability  Assessment  and  Debris-Flow  Initiation  . 25 

Debris-Flow  Rheology . 27 

Debris-Flow  Countermeasures . 27 

Debris-Flow  Magnitude . 28 

Debris-Flow  Modeling  . 28 

Design  and  Construction  of  Low-Water  Crossings . 29 

Sediment  Sampling  and  Analysis . 29 

Annotated  Bibliography  by  U.S.  Geological  Survey  Researchers . 30 

Figures 

1 .  Map  showing  study  area  for  Texas  Department  of  Transportation  research  project  0-4695  .  2 

2.  Graphs  showing  rainfall  hyetographs,  streamflow  hydrographs,  and  sedimenttransport 

rates  for  recorded  floods .  5 

3-6.  Photographs  showing: 

3.  Flume  with  bed  rocks  used  in  experiments  of  Hofland  (Delft  University  of  Technology, 

The  Netherlands,  written  commun.,  2001) .  9 

4.  Side  view  of  flume  and  laser  doppler  used  in  experiments  of  Hofland  (Delft  University 

of  Technology,  The  Netherlands,  written  commun.,  2001 ) .  9 

5.  Experimental  flume  used  for  study  of  pressure  fluctuations  on  and  in  a  subsurface 

gravel  layer .  9 

6.  Different  devices  used  for  study  of  pressure  fluctuations  on  and  in  a  subsurface 

gravel  layer .  10 


IV 


7.  Photograph  and  diagram  showing  riprap  on  hillside  in  Pacifica,  California,  and  riprap 

design  recommendation . 12 

8.  Photographs  showing  stand-alone  flow  breakers  and  flow  breakers  in  debris  dam . 12 

9.  Photograph  and  diagram  showing  debris-flow  dam  with  a  screen  of  steel  beams  and 

open  debris  basin  design  . 13 

10.  Diagram  and  photograph  showing  closed  debris  basin  schematic  and  closed  debris  basin 

on  Harvey  Creek  in  British  Columbia  . 13 

11.  Diagram  and  photograph  showing  deflection  berm  schematic  and  deflection  berm  on 

Boulder  Creek  in  British  Columbia . 14 

12.  Photograph  and  diagram  showing  gabion  deflection  wall,  Tiburon,  California,  and  design 

gabion  deflection  wall . 14 

13.  Diagram  and  photograph  showing  terminal  berm  schematic  and  terminal  berm  on  Cypre 

Creek  in  British  Columbia  . 15 

14.  Photograph  and  diagram  showing  U.S.  Geological  Survey  flume  on  Blue  River  at  Blue  Forest, 
Oregon,  and  schematic  of  the  experiment  to  studythe  response  of  flexible  wire  rope 

barriers  to  debris-flow  charges  . 16 

15.  Photographs  showing  flexible  wire  rope  barriers  with  netting  and  interlocking  rings . 16 

16.  Photographs  showing  gravitational  flow  breaker,  five-layer  grid  flow  breaker,  step-grid 

flow  breaker,  and  column  flow  breaker . 17 

17.  Diagram  showing  flume  test  apparatus  . 17 

18.  Diagram  showing  column  flow  breaker,  beam  flow  breaker,  and  grid  flow  breaker . 18 

19.  Image  showing  processed  Landsat  Thematic  Mapping  scene  for  part  of  western  Arizona  . 21 

Tables 

1.  Grain-size  range  for  Oak  Creek  in  Oregon . 3 

2.  Pavement  and  subpavement  data  for  several  gravel-bed  rivers . 4 

3.  Typical  flow  conditions  in  several  gravel-bed  rivers  . 4 

4.  Hydrologic,  hydraulic,  and  sediment  data  of  recorded  floods . 6 

5.  Debris-flow  management  system  diagram  . 11 


V 


Conversion  Factors 

The  literature  reviewed  uses  inch/pound  (U.S.  customary)  and  SI  (metric)  units  of  measurement. 
The  authors  have  explicitly  chosen  to  remain  consistent  with  the  original  references  in  respect 
to  the  principal  unit  of  measurement  used.  Furthermore,  units  are  not  abbreviated  within  the  text. 

Factors  for  converting  inch/pound  units  to  SI  units: 


Multiply 

By 

To  obtain 

cubic  foot 

0.01832 

cubic  meter 

cubic  foot  per  second 

0.01832 

cubic  meter  per  second 

foot 

0.3048 

meter 

foot  per  second 

0.3048 

meter  per  second 

inch 

25.40 

millimeter 

mile 

1.609 

kilometer 

pound  (mass)  per  foot 

1.488 

kilogram  per  meter 

Temperature  in  degrees  Fahrenheit  (°F)  is  converted  to  degrees  Celsius  (°C)  as  follows: 
°C  =  (°F  -  32)/1 .8 


Factors  for  converting  SI  units  to  inch/pound  units: 


Multiply 

By 

To  obtain 

centimeter 

0.3937 

inch 

cubic  meter 

35.31 

cubic  foot 

cubic  meter  per  second 

35.31 

cubic  foot  per  second 

gram  per  liter 

0.062427 

pound  per  cubic  foot 

kilometer 

0.6214 

mile 

meter 

3.281 

foot 

millimeter 

0.03937 

inch 

millimeter  per  hour 

0.03937 

inch  per  hour 

newton  per  square  meter 

0.020885 

pound  (force  )  per  square  foot 

square  kilometer 

0.3861 

square  mile 

Temperature  in  degrees  Celsius  (°C)  is  converted  to  degrees  Fahrenheit  (°F)  as  follows: 
°F  =  (1.8  x  °C)  +  32 


Blank  Page 


Literature  Review  for  Texas  Department  of 
Transportation  Research  Project  0-4695: 

Guidance  for  Design  in  Areas  of  Extreme 
Bed-Load  Mobility,  Edwards  Plateau,  Texas 

By  Franklin  T.  Heitmuller1,  William  H.  Asquith1,  Xing  Fang2,  David  B.  Thompson3,  and  Keh-Han  Wang4 


Abstract  Introduction 


A  review  of  the  literature  addressing  sediment  transport  in 
gravel-bed  river  systems  and  structures  designed  to  control  bed- 
load  mobility  is  provided  as  part  of  Texas  Department  of  Trans¬ 
portation  research  project  0-4695:  Guidance  for  Design  in 
Areas  of  Extreme  Bed-Load  Mobility.  The  study  area  com¬ 
prises  Edwards,  Kimble,  and  Real  Counties  in  the  Edwards  Pla¬ 
teau  in  central  Texas.  The  primary  focus  of  the  literature  review 
is  on  journal  articles,  edited  volumes,  and  government  publica¬ 
tions.  The  literature  review  provides  an  outline  and  foundation 
for  the  research  project  to  characterize  extreme  bed-load  mobil¬ 
ity  in  rivers  and  streams  across  the  study  area.  The  literature 
review  also  provides  a  basis  for  potential  modifications  to  low- 
water  stream-crossing  design  in  the  study  area.  Major  themes 
within  the  body  of  literature  include  deterministic  sediment 
transport  theory  and  equations,  development  of  methods  to 
measure  and  analyze  fluvial  sediment,  applications  and  devel¬ 
opment  of  theory  in  natural  channels  and  flume  experiments, 
and  recommendations  for  river  management  and  structural 
design.  A  variety  of  methods  are  used  to  study  gravel-bed  river 
characteristics,  including  surveying  channel  geometry,  particle 
tracing,  bed-material  and  bed-load  sampling,  and  developing 
sediment  discharge  rating  curves,  among  others.  Additionally, 
efforts  to  manipulate  the  transport  of  sediment  in  gravel-bed 
rivers  include  construction  of  debris  traps,  velocity  control 
structures,  or  flow  diversion  structures.  Extensive  collection 
and  analysis  of  field  and  geographic  data,  and  subsequent 
parameterization  of  numerical  and  physical  models,  is  the  most 
appropriate  approach  to  understanding  bed-load  transport  and 
structural  design  in  the  Edwards  Plateau. 

1  U.S.  Geological  Survey,  Austin,  Texas 

~  Lamar  University,  Beaumont,  Texas 

3  Texas  Tech  University,  Lubbock,  Texas 

4  University  of  Houston,  Houston,  Texas 


Gravel-bed  rivers  are  inherently  dynamic  features  and  pro¬ 
vide  a  unique  challenge  for  engineers  and  river  management 
specialists.  Activities  in  gravel-bed  rivers  can  span  a  wide  spec¬ 
trum  of  scales,  from  culvert  design  to  river  restoration.  A  suc¬ 
cessful  investigation  of  gravel-bed  rivers  must  consider  flow 
processes,  sediment  characteristics  and  transport,  and  channel 
adjustment  potential.  Extensive  research  of  the  physical  charac¬ 
teristics  and  processes  of  gravel-bed  rivers  has  been  conducted 
by  hydraulic  engineers,  geomorphologists,  hydrologists,  and 
stream  ecologists.  This  research  has  produced  a  large  body  of 
literature  in  numerous  books  and  journals  specific  to  particular 
disciplines.  In  spite  of  the  difficulties  encountered  in  under¬ 
standing  and  predicting  form  and  process  in  gravel-bed  rivers, 
an  appreciable  knowledge  of  the  system  components  is  avail¬ 
able.  In  some  cases,  desired  results  of  engineering  and  river 
management  activities  along  gravel-bed  rivers  have  been 
achieved.  The  successful  implementation  of  a  project  involving 
gravel-bed  rivers  requires  knowledge  of  the  existing  research 
literature. 

A  consortium  of  researchers  at  Texas  Tech  University, 
Lamar  University,  the  University  of  Houston,  and  the  U.S.  Geo¬ 
logical  Survey  (USGS)  was  selected  to  research  the  bed-load 
phenomena  as  part  of  Texas  Department  of  Transportation 
(TxDOT)  research  project  0-4695:  Guidance  for  Design  in 
Areas  of  Extreme  Bed-Load  Mobility.  This  project  was  initiated 
to  address  the  considerable  reconfiguration  of  gravel  during 
high-magnitude  flows  and  the  associated  structural  problems 
and  maintenance  along  low-water  (road)  crossings  (LWCs)  in 
the  study  area.  The  study  area  (fig.  1)  for  the  project  comprises 
Edwards,  Kimble,  and  Real  Counties  in  the  Edwards  Plateau  in 
central  Texas.  A  major  project  objective  was  the  assessment  of 
the  literature.  A  summary  of  the  literature  review  is  provided  in 
this  report.  The  primary  focus  for  the  university  researchers  are 
journal  articles  and  books;  whereas  the  USGS  researchers 
focused  on  government  documents. 
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Figure  1.  Study  area  for  Texas  Department  of  Transportation  research  project  0-4695. 
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Low-water  crossings  are  low-height  structures  that  pro¬ 
vide  acceptable  passage  as  stream  crossings  over  seasonal  riv¬ 
ers  or  streams  with  relatively  low  normal-depth  flow  and  are 
designed  to  allow  flooding  during  high  flow  events.  LWCs  are 
a  common  alternative  to  traditional  bridges  in  rural  areas  and 
along  road  segments  expected  to  receive  low  maintenance 
because  they  generally  are  less  expensive  to  design  and  build 
(Molauskas,  1988;  Pienaar  and  Visser,  1995;  and  Lohnes  and 
others,  2001).  In  areas  that  experience  substantial  bed  mobility, 
LWCs  built  with  no  provision  of  debris-flow  management  dur¬ 
ing  high  streamflow  rates  are  subject  to  structural  erosion,  sed¬ 
iment  clogging,  burial  under  several  feet  of  debris,  and  com¬ 
plete  structural  washout. 

Numerous  LWCs  have  been  built  within  the  Edwards  Pla¬ 
teau  region  of  central  Texas.  During  flood  flow,  LWCs  often  are 
buried  under  several  to  tens  of  feet  of  sediment,  and  the  culverts 
are  clogged  with  gravel-  and  cobble-size  material.  Additional 
damages  include  severe  and  rapid  abrasion  of  the  concrete  cul¬ 
vert  parapet  and  the  pavement.  Structural  failure  on  relatively 
small  LWCs,  which  originally  was  anticipated  to  have  limited 
hydraulic  effect  on  streams,  has  occurred  as  recently  as  the  June 
2004  flood  on  the  Plateau. 

The  purpose  of  this  report  is  to  provide  a  review  of  the  rel¬ 
evant  literature  through  an  annotated  bibliography  and  addi¬ 
tional  description  for  some  of  the  more  informative  citations. 
Furthermore,  this  literature  review  documents  and  synthesizes 
the  extensive  body  of  literature  associated  with  the  characteris¬ 
tics,  entrainment,  and  transport  of  sediment  in  gravel-bed  rivers 
as  well  as  structural  modifications  to  control  sediment  mobility. 
This  work  serves  as  a  reference  for  researchers  and  engineers  to 
increase  their  understanding  of  gravel-bed  river  morphology 
and  sediment  transport  phenomena.  Although  very  few  publica¬ 
tions  directly  address  the  defined  problem  at  LWCs,  many  are 
of  benefit  to  this  investigation,  such  as  structural  designs  to 
minimize  adverse  affects  of  debris  flows.  This  literature  review 
incorporates  different  subjects  applicable  to  the  problem.  Most 
of  the  literature  consists  of  journal  articles,  chapters  in  edited 
volumes,  and  government  documents.  For  readers  new  to  the 
topic,  textbooks  covering  fluvial  geomorphology,  surface-water 
hydrology,  and  the  principles  of  open-channel  flow  also  are  rec¬ 
ommended.  Not  all  references  listed  in  this  report  were 
reviewed  because  of  the  large  volume  of  available  literature. 
References  deemed  most  pertinent  to  the  project  scope  were 
thoroughly  reviewed  and  are  summarized  in  pages  preceding 
the  bibliography. 

Detailed  Literature  Review  of  Selected 
References 

Although  a  wide  variety  of  journals  contain  articles  on  sed¬ 
iment  transport  in  gravel-bed  rivers,  a  handful  distinguish  them¬ 
selves  as  the  preeminent  sources  for  this  type  of  information, 
including  Bulletin  of  the  Geological  Society  of  America',  Earth 
Surface  Processes  and  Landforms',  Geomorphology;  Journal  of 


Geology;  Journal  of  Hydraulic  Engineering  ( Journal  of  the 
Hydraulics  Division );  Journal  of  Hydrology;  Sedimentology; 
and  Water  Resources  Research.  The  emphasis  of  journal  arti¬ 
cles  tends  toward  the  development  of  theories  defining  the  con¬ 
trols  of  gravel  entrainment  and  movement,  resultant  deposi- 
tional  forms,  channel  adjustment  to  the  sediment  regime,  and 
structural  designs  to  control  sediment  mobility.  Five  edited  vol¬ 
umes  on  gravel-bed  rivers  include  Hey  and  others  (1982), 
Thorne  and  others  (1987),  Billi  and  others  (1992),  Klingeman 
and  others  (1998),  and  Mosley  (2001).  These  volumes  contain 
a  variety  of  chapters  discussing  sediment  transport  theory  and 
research  methods.  Finally,  U.S.  government  documents  that 
pertain  to  sediment  transport  in  gravel-bed  rivers  are  most  fre¬ 
quently  published  by  the  USGS  and  the  U.S.  Army  Corps  of 
Engineers.  Early  government  publications  introduced  theories 
of  bed- load  transport  and  channel  adjustment,  and  more  recent 
publications  discuss  project  applications  involving  bed-load 
transport  and  provide  methodological  guidance. 

Literature  Reviews  by  Lamar  University  Researchers 

SedimentTransport  Mechanism  in  Gravel-Bed  Rivers 

The  definition  of  “sediment  transport”  is  the  transport  of 
granular  particles  by  fluid.  As  such,  it  embodies  a  type  of  two- 
phase  flow,  where  one  is  fluid  and  the  other  is  a  solid  composed 
of  granular  particles.  Flow  condition,  grain-size  distribution, 
and  boundary  shear  stress  are  the  main  factors  affecting  the  sed¬ 
iment  transport  mechanism.  Gravel-bed  streams  usually  have  a 
bed  surface,  referred  to  as  “pavement,”  markedly  coarser  than 
the  substrate  (Milhous,  1973;  Parker,  Klingeman,  and  McLean, 
1982).  Although  the  name  may  be  misleading,  pavement  is 
mobile  and  differs  from  an  immobile  armored  layer.  The  grain- 
size  distribution  for  Oak  Creek  in  Oregon  is  listed  in  table  1 ,  and 


Table  1.  G  rain-size  range  for  Oak  Creek  in  Oregon  (from  Parker, 
Dhamotharan,  and  Stefan,  1982). 

[mm,  millimeters;  fj,  subpavement;  fpj,  surface  pavement;  <,  less  than;  — ,  no 
data] 


Range 

(mm) 

Grain  diameter 
(mm) 

fi 

(percent) 

fpi 

(percent) 

203-102 

152 

0 

7.7 

102-76 

88.9 

4 

18.4 

76-51 

63.5 

15 

26.4 

51-38 

44.4 

10 

13.0 

38-25 

31.8 

15 

12.8 

25-19 

22.2 

6.5 

6.7 

19-9.5 

14.3 

13 

6.1 

9.5^1.8 

7.14 

11.5 

3.5 

4. 8-2.4 

3.57 

10 

1.8 

2.4-1. 2 

1.79 

8 

1.1 

1.2-.59 

.890 

4 

.9 

.59-30 

.444 

1.9 

1.0 

.30-.  15 

.223 

.7 

.3 

.15-.  74 

.112 

.2 

.2 

<.74 

- 

.2 

.1 
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the  data  confirm  that  surface  pavement  has  much  coarser 
sediment  than  subpavement.  Table  2  lists  pavement  and  sub¬ 
pavement  data  for  several  gravel-bed  rivers;  median  grain 
diameters  of  pavement  ranged  from  44  to  76  millimeters. 
Typical  flow  conditions  in  several  gravel-bed  rivers  previously 
examined  are  listed  in  table  3  (Hollingshead,  1971;  Jones  and 
Seitz,  1980).  The  beds  of  paved  streams  require  only  infrequent 
flow  conditions  to  activate  the  bed  motion;  therefore,  the 
bed  is  active  for  at  least  several  days  in  most  years  (Parker, 
Klingeman,  and  McLean,  1982).  Milhous  (1973)  determined 
that  the  pavement  is  first  “broken”  at  a  critical  water  discharge 
of  about  1 . 1  cubic  meters  per  second  (about  40  cubic  feet 
per  second).  During  such  periods,  all  or  nearly  all  sizes,  includ¬ 
ing  the  coarsest  pavement  grains,  often  are  found  in  motion 
(Hollingshead,  1971;  Milhous,  1973).  Many  gravel-bed  streams 
can  transport  essentially  all  available  grain  sizes  during  suffi¬ 
ciently  large  floods  and  yet  show  the  same  pavement  at  low 
flow  (Parker  and  Klingeman,  1982).  Parker  and  Klingeman 
(1982)  propose  two  hypotheses: 


Table  2.  Pavement  and  subpavement  data  for  several  gravel-bed 
rivers  (from  Parker,  Dhamotharan,  and  Stefan,  1982). 


[mm,  millimeters;  D90,  grain  size  at  which  90  percent  of  sediment  is  finer;  D50, 
grain  size  at  which  50  percent  of  sediment  is  finer;  *,  sand  bars  on  top  of  the 
pavement  are  frequent.] 


Stream 

name 

Pave¬ 

ment 

D90 

(mm) 

Pave¬ 

ment 

D50 

(mm) 

Bed 

material 

D50 

(mm) 

Percent 
sand 
in  bed 

Oak  Creek 

86 

54 

20 

12 

Elbow  River 

132 

76 

28 

67 

Snake  River 

137 

54 

27 

14* 

Clearwater  River 

143 

72 

18 

19* 

Vedder  River 

90 

44 

19 

16 

Table  3.  Typical  flow  conditions  in  several  gravel-bed  rivers 
(from  Parker,  Dhamotharan,  and  Stefan,  1982). 

[m.  meters;  ft3/s,  cubic  feet  per  second] 


Stream 

name 

Water- 

surface 

slope 

Water- 

surface 

width 

(m) 

Average 
depth  or 
hydraulic 
radius 
(m) 

Discharge 

(ft3/s) 

Oak  Creek 

0.0097-0.0108 

5.0-6. 1 

0.31-0.45 

1.16-3.40 

Elbow  River 

.00745 

38-49 

0.64-0.86 

35-109 

Snake  River 

.00094-00121 

184-198 

4. 6-5.9 

2,200-3,500 

Clearwater 

River 

.00035-00062 

142-149 

4. 9-6.4 

1,500-3,100 

Vedder 

River 

.00195 

85-90 

1.34-1.66 

216-370 

1.  During  floods  capable  of  moving  pavement  grains,  the 
pavement  ceases  to  exist  or  is  greatly  modified.  As  flow  wanes, 
it  reforms. 

2.  Pavement  is  present  even  during  floods  capable  of  mov¬ 
ing  all  available  grains.  Its  structure  during  such  floods  differs 
little  from  that  at  low  flow. 

Parker,  Klingeman,  and  McLean  (1982)  conducted  exper¬ 
iments  that  show  that  the  second  hypothesis  is  the  correct 
one.  Pavement  formed  readily  and  could  be  maintained  indefi¬ 
nitely  under  constant  flows,  even  though  all  grain  sizes  were 
represented  in  the  bed  load.  In  other  words,  pavement  is  a 
mobile  bed  phenomena  (Parker  and  Klingeman,  1982).  Parker 
and  Klingeman  (1982)  state  that  pavement  is  a  regulator  that 
enables  a  stream  to  transport  the  coarse  one-half  and  fine  one- 
half  of  its  bed-load  supply  at  equal  rate.  More  specially,  pave¬ 
ment  forms  just  so  as  to  render  all  available  grain  sizes  of  nearly 
equal  mobility.  Parker  and  Klingeman  (1982)  state  that  if  the 
size  distributions  of  bed  load  and  bed  material  are  almost  simi¬ 
lar,  the  coarse  one-half  of  the  subpavement  could  move  through 
the  reach  at  the  same  rate  as  the  fine  one-half  if  a  reach  is  to  be 
in  equilibrium. 

Parker  and  Klingeman  (1982)  consider  a  phenomenon 
called  “hiding,”  previously  noted  by  Einstein  (1950)  and 
Egiazaroff  (1965).  The  “microscopic”  hiding  depends  on  the 
relative  placement  of  individual  grains.  The  isolated  protrusions 
offered  by  coarser  surface  grains  in  a  mixture  render  them  more 
mobile  than  they  are  in  a  uniform  sediment  composed  exclu¬ 
sively  of  the  coarse  size  (Parker  and  Klingeman,  1982).  It 
means  that  in  a  mixture,  due  to  microscopic  hiding,  finer  surface 
grains  are  less  mobile,  and  coarse  grains  are  more  mobile  in 
comparison  to  a  uniform  sediment  condition.  Microscopic  hid¬ 
ing  helps  reduce  the  intrinsic  difference  in  mobility  between 
coarser  and  finer  surface  grains  (Parker  and  Klingeman,  1982). 
The  surface  pavement  also  provides  a  macroscopic  hiding 
effect,  so  that  gravel-bed  streams  eliminate  the  remaining 
mobility  difference  by  means  of  the  pavement  itself.  No  matter 
what  percentage  of  finer  grains  of  a  given  size  are  stored  in  the 
subpavement,  none  can  be  entrained  if  none  are  contained  in  the 
pavement  (Parker  and  Klingeman,  1982).  The  authors  also 
explain  the  vertical  winnowing  of  fine  particles.  Each  time  a 
large  grain  is  dislodged,  it  leaves  a  “hole”  of  comparable  size  in 
the  bed.  Small  grains  that  fall  into  the  hole  might  work  their  way 
below  the  pavement  and  reduce  their  probability  of  re-erosion. 
This  was  also  described  by  Milhous  (1973). 

Parker,  Klingeman,  and  McLean  (1982)  use  similarity 
analysis  to  study  the  size  distribution  in  paved  gravel-bed 
streams  and  conclude  that  all  grain-size  ranges  are  of  about 
equal  transportability  when  the  critical  condition  for  breaking 
the  pavement  is  exceeded.  When  discharge  exceeds  the  critical 
value,  more  and  more  of  the  pavement  is  moved,  exposing  the 
finer  gravels  below  and  making  them  available  for  transport. 
Bed-load  transport  is  governed  by  hydraulic  conditions  rather 
than  availability. 

Parker,  Dhamotharan,  and  Stefan  (1982)  performed 
laboratory  experiments  with  a  poorly  sorted  gravel  bed  and 
intended  to  provide  a  mechanistic  model  for  two  paved  streams 
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Figure  2.  Rainfall  hyetographs,  streamflow  hydrographs,  and  sodimont  transport  rates  for  recorded  floods  (from  Billi  and  others,  1998). 


for  which  bed  load  had  been  measured.  These  streams  are  the 
Elbow  River  in  Alberta,  Canada  (Hollingshead,  1968)  and  Oak 
Creek  in  Oregon  (Milhous,  1973).  Both  have  a  steep  down¬ 
stream  water-surface  slope  (S)  (0.0075,  Elbow  River;  0.010, 
Oak  Creek).  The  results  reveal  that  pavement  is  a  mobile  bed 
phenomenon,  but  the  pavement  can  coexist  with  the  motion  of 
all  available  grain  sizes.  Because  of  its  episodic  or  sporadic 
motion,  at  any  given  time  only  a  small  percentage  of  surface 
grains  are  actually  in  motion.  Therefore,  particles  in  the  imme¬ 
diate  subpavement  are  only  occasionally  mobilized  and  deeper 
particles  are  rarely  moved,  even  during  floods.  It  was  found  that 
pavement  imbrication  is  ubiquitous,  even  though  pavement  and 
bed  load  exchange  with  each  other;  at  higher  flow  there  was 
weak  formation  of  bars. 

Kuhnle  (1992)  documents  bed-load  transport  during  ris¬ 
ing  and  falling  stages  on  two  small  streams.  At  higher  flow 
strengths,  mean  bed-load  transport  rates  were  greater  during  ris¬ 
ing  stages  than  during  falling  stages.  One  of  the  streams  showed 
evidence  for  greater  transport  rates  for  low  flows  as  the  stage 
declined.  However,  experiments  from  Sutherland  and  cowork¬ 
ers  (Griffiths  and  Sutherland,  1977;  Bell  and  Sutherland,  1983; 
and  Phillips  and  Sutherland,  1989,  1990)  show  that  unsteady 
flow  does  not  affect  the  rate  of  bed-load  transport  unless  the  bed 
is  degrading. 


Billi  and  others  (1998)  investigate  channel  morphology 
processes  and  bed  load  in  a  high  altitude  alpine  torrent  in  the 
Rio  Cordon  catchment.  The  authors  use  sediment  transport 
measurements  to  study  relations  between  slope  erosion  pro¬ 
cesses,  bed  morphology  adjustment,  and  sediment  transport. 
The  study  area  has  the  following  characteristics:  watershed 
area,  52.0  square  kilometers;  main  channel  slope,  17  percent; 
annual  precipitation,  1,100  millimeters;  and  average  annual 
temperature,  2.0  °C. 

Table  4  lists  sediment  transport  values  from  eight  floods 
(Billi  and  others,  1998).  Peak  discharges  range  from  1.8  to  10.4 
cubic  meters  per  second  (or  64  to  370  cubic  feet  per  second), 
two  critical  discharges  for  bed-load  entrainment  occur  at  2 
cubic  meters  per  second  (or  70  cubic  feet  per  second),  and 
median  particle  size  of  coarse  bed  load  (D50)  ranges  from  25.8 
to  82.6  millimeters  (1.0  to  3.3  inches).  Figure  2  depicts  hyeto¬ 
graphs,  hydrographs,  coarse  bed  load,  and  suspended  load  dur¬ 
ing  the  floods.  The  flood  on  September  14,  1994,  the  most 
severe  flood  recorded,  was  caused  mostly  by  high  intensity  rain¬ 
fall  (maximum  rainfall  rate  7.2  millimeters/5  minutes,  16.4  mil- 
limeters/15  minutes,  and  25.3  millimeters/30  minutes).  The 
recurrence  interval  was  assessed  to  be  between  35  and  50  years. 
Two  flood  waves  occurred,  with  the  second  having  a  higher 
peak  discharge  (10  cubic  meters  per  second).  The  volume  of 
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Table  4.  Hydrologic,  hydraulic,  and  sediment  data  of  recorded  floods  (modified  from  Billi  and  others,  1998). 


[  m 3 ,  cubic  meters;  m3/s,  cubic  meters  per  second;  h,  hours;  m3/h,  cubic  meters  per  hour;  mm,  millimeters;  g/L,  grams  per  liter;  (1)(2),  first  and  second  part  of 
flood;  — ,  not  available;  [  ],  assessed  value] 


Sediment 

transport 

characteristic1 

Oct.  11, 
1987 

July  3, 
1989 

June  17, 
1991 

Oct.  5, 
1992 

Oct.  2, 
1993 

May  19, 
1994 

July  18, 
1994 

Sept.  14, 
1994 
(1) 

Sept.  14, 
1994 
(2) 

Re(103m3) 

79.90 

103.36 

57.89 

21.52 

30.69 

5.41 

0 

4.71 

21.93 

Qp(m3/s) 

5.15 

4.4 

4.0 

2.9 

4.3 

1.80 

1.8 

3.74 

10.4 

Yb(m3) 

50.0 

85.0 

39.0 

9.3 

10.2 

1.00 

0 

10 

890 

T(h) 

8 

27 

20 

10 

[6] 

12 

0 

1.42 

2.75 

Qb(m3/h) 

6.25 

3.15 

1.95 

.93 

[1.7] 

.08 

0 

7 

323.6 

D5o(mm) 

- 

82.5 

35.5 

31.1 

43.1 

25.8 

- 

- 

82.6 

D84(mm) 

- 

152.0 

64.0 

77.0 

96.0 

37.0 

- 

- 

168.0 

D90(mm) 

- 

173.0 

84.0 

117.0 

133.0 

43.5 

- 

- 

207.0 

Qcri(m3/s) 

1.8 

2.2 

[2.0] 

1.9 

[2.3] 

1.60 

- 

1.8 

1.8 

Qcr2(m3/s) 

3.8 

2.7 

[2.4] 

2.1 

[3.7] 

1.61 

- 

1.8 

3.3 

Cmax(g/L) 

- 

- 

- 

1.2 

1.65 

.1 

7.8 

30 

58 

^Definition  of  symbols: 


Re  =  effective  runoff; 

Qp  =  peak  discharge; 

Yb  =  deposited  volume  of  coarse  bed  load; 

T  =  bed-load  transport  duration; 

Qb  =  average  transport  rate  of  coarse  bed  load; 
D50  =  median  particle  size  of  coarse  bed  load; 


D84>  D90 

Qcrl 

Qcr2 

C 

'-'max 


=  particle  sizes  for  which  84  and  90  percent,  respectively,  of  coarse 
bed-load  material  is  finer; 

=  discharge  at  bed-load  entrainment; 

=  discharge  at  bed-load  transport  end; 

=  maximum  suspended  sediment  concentration. 


coarse  sediment  transported  was  about  900  cubic  meters,  much 
larger  than  the  volumes  measured  for  the  previous  floods,  which 
commonly  were  less  than  100  cubic  meters.  It  also  was 
observed  that  bed-load  transport  occurs  at  two  different  scales 
of  magnitude.  Typical  bed  loads  are  associated  with  peak  flow 
discharges  less  than  5  cubic  meters  per  second  with  a  return 
period  not  exceeding  5  years  and  bed-load  transport  rates  of 
7  cubic  meters  per  hour.  During  larger  floods  with  return  period 
exceeding  30  years,  the  bed-load  transport  rate  is  as  much  as 
324  cubic  meters  per  hour,  or  50  times  the  transport  of  a  typical 
bed  load. 

Sechet  and  Le  Guennec  (1999)  introduce  the  bursting 
phenomenon  concept  associated  with  the  incipient  motion  of 
the  solid  particles  in  bed-load  transport.  They  conducted  an 
experiment  to  visualize  the  displacement  of  solid  particles  on 
the  smooth  bottom  of  a  hydraulic  channel  and  calculated  trajec¬ 
tories  with  the  laser  doppler  anemometry  measurement  of  the 
instantaneous  velocity  near  the  bottom.  Sechet  and  Le  Guennec 
(1999)  conclude  that  there  coexist  two  transport  modes  at  the 
wall:  one  is  transport  by  sweeps  and  the  other  is  through  ejec¬ 
tion.  Transport  by  the  ejection  has  the  dominating  effect,  and 
the  particle  transport  by  sweeps  needs  further  investigation. 

McEwan  and  Heald  (2001)  examined  the  stability  of  ran¬ 
domly  deposited  sediment  beds  by  using  a  discrete  particle 
model  in  which  individual  grains  are  represented  by  spheres. 
The  results  indicate  that  the  threshold  or  critical  entrainment 
shear  stress  for  cohesionless  uniform  size  sediments  in  flat  beds 
cannot  be  adequately  described  by  a  single-value  parameter; 


rather,  it  is  best  represented  by  a  distribution  of  values.  They 
state  that  a  distinct  threshold  for  the  initiation  of  motion  does 
not  exist  but  a  value  of  0.06  traditionally  has  represented  the 
threshold  condition  for  the  uniform  size  sediment.  There  are 
two  types  of  sheltering  that  occur  in  the  sediment  bed:  direct 
sheltering  and  remote  sheltering.  In  the  case  of  direct  sheltering, 
particles  upon  which  a  particular  grain  is  resting  reduce  the 
grain’s  effective  area  and,  hence,  alter  the  fluid  force  acting 
upon  it.  This  condition  will  remain  as  long  as  the  grain  is  static 
in  its  resting  condition.  Remote  sheltering  is  a  consequence  of 
non-uniformity  of  the  bed  flow  that  occurs  over  complex  rough¬ 
ness.  For  remote  sheltering,  other  upstream  particles,  poten¬ 
tially  located  a  number  of  diameters  upstream,  modify  the  fluid 
velocity  in  the  vicinity  of  the  grain  and,  hence,  influence  the 
fluid  forces  acting  on  it.  They  also  state  that  remote  sheltering 
and  direct  sheltering,  the  supporting  grain  arrangement,  and  the 
influence  of  remote  sheltering  by  upstream  grains  dominate  the 
critical  shear  stress  distribution.  Direct  sheltering  is  less  influ¬ 
ential,  comparatively. 

Parker  and  Toro-Escobar  (2002)  conclude  that  there 
are  two  hypotheses  on  the  concept  of  equal  mobility  of  gravel 
transport.  The  weak  hypothesis  states  that  to  transport  the 
coarse  one-half  of  the  mean  annual  gravel  load  through  a 
river  reach  at  a  rate  equal  to  the  fine  one-half,  the  coarse 
material  must  be  overrepresented  on  the  bed  surface,  giving 
rise  to  mobile  bed  armor.  The  strong  hypothesis  states  that  the 
grain- size  distribution  of  the  gravel  portion  of  the  bed  load 
should  be  similar  to  that  of  the  substrate  and  finer  than  that  of 
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the  surface  layer.  A  rigorous  test  in  the  laboratory  under  some 
controlled  conditions,  which  is  almost  similar  to  the  field  con¬ 
dition,  was  conducted.  The  results  of  the  experiment  strongly 
support  the  weak  hypothesis;  the  surface  layer  was  substantially 
coarser  than  the  bed-load  material.  The  strong  hypothesis  is  also 
supported  by  the  experiment,  in  that  the  size  distribution  of  the 
“gravel”  substrate  was  nearly  identical  to  the  size  distribution  of 
the  “gravel”  portion  of  the  load  (Parker  and  Toro-Escobar, 
2002).  The  gravel  substrate  was,  however,  slightly  but  measur¬ 
ably  coarser  than  the  gravel  load,  a  tendency  that  has  been 
observed  in  the  field  (Parker  and  Toro-Escobar,  2002). 

Hydraulic  Geometry  of  Gravel-Bed  Rivers 

Parker  (1979)  discusses  active  gravel  rivers  that  satisfy 
the  following  criteria:  (1)  channel  is  straight  and  laterally  sym¬ 
metrical;  (2)  perimeter  is  composed  of  loose,  coarse  gravel  that 
is  of  similar  composition  in  both  bed  and  banks  and  that  cannot 
be  suspended;  (3)  most  particles,  including  those  much  coarser 
than  median  size,  are  mobile  at  bankfull  or  dominant  condi¬ 
tions;  and  (4)  self-formed  channel  has  stable  width  and  is  in 
grade.  Parker  (1979)  reveals  that  wide  gravel-bed  rivers  have 
larger  loads  in  comparison  to  the  narrow  channels.  The  relations 
regarding  the  effect  of  bed-load  material  on  channel  geometry 
and  vice  versa  are  given  by  the  following  equations. 


transport  rate  to  hydraulic  and  sedimentological  quantities  from 
field  data,  laboratory  data,  or  mechanical/physical  principals. 
Many  equations  have  been  modified  in  light  of  their  perfor¬ 
mance  against  data  not  included  in  the  initial  analysis.  Most  of 
the  equations  have  been  assessed  tentatively  on  the  basis  of  lim¬ 
ited  amount  of  field  data.  No  completely  adequate  methodology 
has  been  developed  to  assist  the  investigator  in  selecting  the 
best  equation  for  a  particular  hydraulic  and  sediment  condition. 
If  no  sediment  transport  information  is  available  for  a  river, 
none  of  the  selected  equations  are  capable  of  generally  predict¬ 
ing  bed-load  transport  in  gravel-bed  rivers  (Gomez  and  Church, 
1989).  Several  recent  equations  that  have  been  reviewed  and 
summarized  follow. 

Parker,  Klingeman,  and  McLean  (1982)  analyzed  data 
from  Oak  Creek  and  conclude  that  for  conditions  in  excess  of 
the  critical  stress  of  the  pavement,  the  bed-load  size  distribution 
is  to  a  first-order  approximation  the  same  as  the  subpavement 
size  distribution.  Empirical  relations  for  total  bed-load  material 
as  a  function  of  shear  stress  and  subpavement  median  size  were 
developed  for  poorly  sorted  gravel-bed  streams.  The  dimen¬ 
sionless  bed  load  ( W* )  and  shear  stress  (x*  )  for  the  ith  size 
range  are  defined  as: 

yy*  _  \ gClBi  _  ^ ClBi  _  ClBi 
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H*c  =  0.05535  1013 ,  (3) 


T  =  p gdS , 


(7) 


where  Qs  is  the  volumetric  gravel  discharge  of  bed  material; 
A=  ps/(p  -  1 )  is  the  submerged  specific  gravity  of  the  sedi¬ 
ment  with  p  and  p  as  densities  of  water  and  sediment;  g  is  the 
acceleration  of  gravity;  Ds0  is  the  median  grain  diameter  at  the 

surface  pavement  layer;  H*c=  II c/ D50  is  the  dimensionless 

* 

depth  at  the  channel  central  region;  B  =  B/D50  is  the  dimen¬ 
sionless  channel  width;  and  Q  =  Q /  ( JXgD 50D"50)  is  the 
dimensionless  water  discharge. 

SedimentTransport  Equations  for  Gravel-Bed  Rivers 

Four  principal  approaches  have  emerged  in  the  design  of 
bed-load  transport  equations  (Gomez  and  Church,  1989),  on  the 
basis  of  bed  shear  stress  (DuBoys,  1879),  stream  discharge 
(Schoklitsch,  1934),  stochastic  functions  for  sediment  move¬ 
ment  (Einstein,  1950),  or  stream  power  (Bagnold,  1980).  A  for¬ 
mal  or  theoretical  relation  was  derived  to  link  the  bed-load 


where  qBi  is  volumetric  bed-load  transport  rate  per  unit  channel 
width;  f\  is  fraction  of  the  bed  material;  D!  is  typical  (mean) 
grain  size,  for  i  size  range  of  bed  load;  d  is  the  hydraulics  radius 
or  cross-sectionally  averaged  depth;  X  =  (ps/p-l)  is  the 
submerged  specific  gravity  of  the  sediment;  and  S  is  the  down¬ 
stream  slope  of  energy  grade  line.  After  integration  of  all  size 
ranges,  the  total  bed-load  relation  is  only  a  function  of  the  sub¬ 
pavement  median  grain  size. 

W*  =  0.0025 exp [  14.2 (c|)50-  1)  -  9.28(<|)50  -  l)2]  for 

0.95  <  <|>50  <  1.65  ,  and  (8) 

W*  =  11.2(1  -  0.822/c|>50)4'5  for  cf)50  >  1.65  .  (9) 

where  (f»50  =  x*50/ x*50  is  the  normalized  Shields  stress;  the 
reference  value  x*50  was  0.0876  (Parker,  Dhamotharan,  and 
Stefan,  1982),  which  is  associated  with  the  reference  W*  of 
0.002.  The  substrate-based  bed-load  relation  above  strictly 
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applies  only  to  near-equilibrium  mobile-bed  conditions.  It 
is  useful  for  rough  estimates  of  gravel  yield  even  in  disequilib¬ 
rium  conditions  (Parker,  1990).  A  substrate-based  relation  is 
not  suited  to  the  computation  of  selective  sorting  driven  by 
disequilibrium. 

Parker  and  Klingeman  (1982)  introduce  the  relative 
mobility  of  different  uniform  materials.  The  relative  mobility 
analysis  reveals  that,  for  the  range  of  stresses  typical  of  field 
gravel-bed  streams,  uniform  coarser  material  is  substantially 
less  mobile  than  uniform  finer  material  at  the  same  fluid  stress. 
However,  when  dimensionless  Shields  stress  is  greater  than  the 
critical  Shields  stress,  the  relative  mobility  is  more  or  less  inde¬ 
pendent  of  the  grain-size  ratio;  therefore,  approximate  equal 
mobility  is  attained  (Parker  and  Klingeman,  1982).  The  dimen¬ 
sionless  bed  load  was  normalized  to  be  a  function  of  the  normal¬ 
ized  Shields  stress  as: 


W*/W*ref  =  G(4>)  =  G(  tV<) 

5.6x103(  1  -  0.853/ (|))4'5 ,  and  (10) 

I,!  =  1.17x*  ,  (II) 

where  x*  is  the  critical  Shields  stress,  and  W*rej-  was  0.002  for 
Oak  Creek  data  (Parker  and  Klingeman,  1982).  The  total  bed¬ 
load  computation  is  given  as: 


qBi  =  VfiJgidsy'^KejGW/^g)  •  (12) 

Parker  (1990)  derives  the  surface-based  bed-load 
transport  equation  for  gravel-bed  rivers  by  using  the  data  from 
Oak  Creek  in  Oregon,  because  mobilization  of  the  grains  takes 
place  from  the  exposed  bed  surface  and  is  caused  by  the  fluid 
force  on  the  exposed  grains.  Substrate  can  participate  as  bed 
load  only  under  the  condition  of  local  or  global  scour.  The  sur¬ 
face-based  bed-load  relation  developed  by  Parker  (1990)  is: 
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G(tf) )  =  (])14'“  for  §  <  1 ,  (18) 

where  W*sj  is  the  dimensionless  bed  load  for  the  ith  size  range 
of  surface  grains;  gs  is  a  “reduced”  hiding  function;  Dsg  is  the 
geometric  mean  surface  size;  and  coo(c|)sgo)  is  a  straining  func¬ 
tion  ranging  from  1.01 1  for  very  coarse  grains  to  0.453,  for 
which  the  equilibrium  surface  layer  is  considerably  finer  and 
poorly  sorted. 

Gomez  and  Church  (1989)  test  the  performance  of  the  12 
bed-load  sediment  transport  equations  for  gravel-bed  channels 
with  seven  sets  of  data  (four  sets  of  river  data  and  three  sets  of 
flume  data).  Gomez  and  Church  (1989)  conclude  that  none  of 
the  equations  performed  consistently  well.  Of  seven  complete 
comparisons,  only  three  equations  account  for  successful  mean 
comparison:  (1)  Schoklitsch  (1934),  (2)  Bagnold  (1980),  and 
(3)  Einstein  (1950).  The  Bagnold  (1980)  formula  is  predomi¬ 
nantly  more  successful  in  comparison  to  others.  To  understand 
the  magnitude  of  bed-load  transport  on  the  basis  of  limited 
hydraulic  information,  the  stream  power  equation  should  be 
used  (Gomez  and  Church,  1989).  When  local  transport  esti¬ 
mates  are  required  and  local  hydraulic  information  is  available, 
an  equation  should  be  selected  that  is  sensitive  to  bed  state  or 
grain-size  distribution,  including  the  Einstein  ( 1950),  Parker 
(1979),  and  Ackers  and  White  (1973)  equations. 

Bakke  and  others  (1999)  calibrate  the  Parker  and 
Klingeman  (1982)  model  to  minimize  the  variance  and  bias. 
The  model  carries  several  assumptions;  one  assumption  is  that 
the  bed  load  is  in  equilibrium  with  the  bed  material,  and  the  sys¬ 
tem  is  neither  aggrading  nor  degrading.  The  pavement  layer  is 
formed  to  establish  an  equilibrium  between  the  bed  load  and  the 
bed  material.  Parker  and  Klingeman  (1982)  use  the  “hiding  fac¬ 
tor”  to  predict  bed  load.  This  helps  increase  the  critical  shear 
stress  required  for  motion  of  the  smaller  particles  compared  to 
that  for  a  uniform  bed,  which  slows  down  the  transfer  of  smaller 
particles.  The  hiding  force  is  described  by  the  simple  relation: 


T/-50 


h 


(19) 


where  x*(  and  x*50  are  the  reference  Shields  stresses  associated 
with  subpavement  Dt  and  Ds0 ,  and  h  =  -0.982  for  Oak  Creek. 
Bakke  and  others  (1999)  indicate  that  h  ranged  from  -0.907 
to  -0.974  on  the  basis  of  sediment-transport  data  from  seven 
Klamath  River  Basin  streams.  The  authors  suggest  a  calibration 
of  h  with  site-specific  data  for  application  of  the  Parker- 
Klingeman  equation.  Lopez  and  Garcia  (2001)  used  a  probabi¬ 
listic  approach  to  analyze  sediment  transport. 


G((f))  =  5474(1  -0.853/c|))4'5  for  <|>  >  1.59,  (16) 

G(<»  =  exp [  14.2(cf)  -  1)  -  9.28(c))  -  l)2]  for 

1<<|><  1.59,  and  (17) 


Laboratory  Experiments  for  Gravel  Transport 

B.  Hofland,  a  Ph.D.  student  in  the  Fluid  Mechanics 
Section,  Faculty  of  Civil  Engineering  and  Geosciences,  Delft 
University  of  Technology,  The  Netherlands,  conducted  a 
series  of  laboratory  experiments  to  study  pressure  and  velocity 
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fluctuations  around  a  granular-bed  element.  In  the  study,  rocks 
are  used  to  prevent  scouring  of  the  underlying  bed  material 
(fig.  3).  The  stability  of  stones  in  such  a  granular  filter  is  usually 
evaluated  by  means  of  the  non-dimensional  bed  shear  stress  or 
Shields  factor.  This  factor  is  a  measure  of  the  mobility  of  the 
bed;  if  higher  than  a  certain  threshold,  the  bed  protection  is  sub¬ 
ject  to  damage  (B.  Hofland,  Delft  University  of  Technology, 
The  Netherlands,  written  commun.,  2001).  However,  if  the  flow 
is  not  an  equilibrium  boundary-layer  flow,  this  approach  (using 
only  the  mean  bed  shear  stress)  does  not  apply.  A  good  example 
is  the  flow  behind  a  backward-facing  step.  At  the  reattachment 
point,  the  mean  bed  shear  stress  vanishes.  Nevertheless,  many 
stones  are  entrained  at  this  location,  because  of  the  increased 
level  of  turbulence  (B.  Hofland,  Delft  University  of  Technol¬ 
ogy,  The  Netherlands,  written  commun.,  2001).  Holland's 
research  was  aimed  at  gaining  more  insight  into  the  interaction 
between  the  turbulence  structure  and  the  forces  on  the  stones 
to  develop  new  (physically  based)  design  criteria  to  select 
appropriate  rock  sizes  to  prevent  scouring  of  the  underlying  bed 
material. 


Figure  4.  Side  view  of  flume  and  laser  doppler  used  in  experi¬ 
ments  of  B.  Hofland  (Delft  University  of  Technology,  The 
Netherlands,  written  commun.,  2001). 


Figure  3.  Flume  with  bed  rocks  used  in  experiments  of  B.  Hofland 
(Delft  University  of  Technology,  The  Netherlands,  written 
commun.,  2001). 

To  get  a  good  picture  of  the  flow  processes  near  and  in  a 
granular  bed,  a  number  of  accurate  pressure  sensors  were  used 
with  a  laser  doppler  (fig.  4)  (B.  Hofland,  Delft  University  of 
Technology,  The  Netherlands,  written  commun.,  2001).  The 
pressure  sensors  give  a  detailed  view  of  the  pressure  field.  In 
combination  with  velocity  measurements,  the  position  and  kind 
of  turbulence  flow  structures  that  create  the  extreme  forces  on 
stones  are  determined. 

Martin  Detert  and  several  researchers  at  the  Institute  for 
Hydromechanics,  University  of  Karlsruhe;  the  Interdisciplinary 
Center  for  Scientific  Computing,  University  of  Heidelberg;  and 
the  Federal  Waterways  Engineering  and  Research  Institute, 


Germany,  conducted  laboratory  experiments  to  study  pressure 
fluctuations  on  and  in  a  subsurface  gravel  layer  caused  by  tur¬ 
bulent  open-channel  flow  (fig.  5).  The  experiments  were  made 
to  investigate  riverbed  stability  and  to  quantify  the  influence  of 
turbulence  in  the  free  surface  flow  to  the  pore  flow.  Pore  flow 
was  measured  by  tracking  tracer  particles  with  two  flexible 
optic  endoscopes  (shown  at  left  in  fig.  6).  Both  velocity  and 
pressure  of  pore  flow  were  measured.  Additional  micropressure 
sensors  allowed  simultaneous  gaging  of  pressure  fluctuations 
within  the  gravel  layer  (shown  at  center  in  fig.  6).  Results  were 
presented  in  the  Conference  of  River  Flow  2004,  at  Naples, 
Italy,  lune  23-25,  2004  (Detert  and  others,  2004;  Jehle  and 
others,  2004). 


Figure  5.  Experimental  flume  for  study  of  pressure  fluctuations 
on  and  in  a  subsurface  gravel  layer  (picture  courtesy  of  Martin 
Detert,  University  of  Karlsruhe). 
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Figure  6.  Different  devices  used  for  study  of  pressure  fluctua¬ 
tions  on  and  in  a  subsurface  gravel  layer. 

Lane  and  others  (2004)  document  both  an  experimental 
and  numerical  study  on  coherent  flow  structures  in  shallow 
gravel-bed  rivers.  The  authors  used  three-dimensional 
computational  fluid  dynamics,  particle  image  velocimetry,  and 
acoustic  Doppler  velocimetry  to  explore  the  nature  of  flow 
structures  in  a  flume  study  of  shallow  gravel-bed  rivers.  This 
research  is  based  on  a  new  algorithm  for  representing  very 
complex  topography  in  structured  grids.  The  model  study 
effectively  represented  the  changing  flow  field  around  an 
individual  bed  particle,  and  it  has  major  potential  benefits 
toward  the  understanding  of  flow-form  interactions  and  the 
associated  implications  for  sediment  transport  and  shallow  flow 
ecohydraulics. 

Literature  Reviews  byTexasTech  University 
Researchers 

Researchers  at  Texas  Tech  University  made  contributions 
to  the  literature  review  by  coordinating  responsibilities  among 
the  other  research  entities,  communicating  with  TxDOT,  and 
making  an  extensive,  but  unsuccessful,  search  for  publications 
directly  addressing  problems  associated  with  sediment  mobility 
at  LWCs.  The  lack  of  literature  addressing  LWCs  could  indicate 
either  that  problems  may  be  unique  to  the  Edwards  Plateau  or 
that  research  has  not  been  conducted  to  understand  the  effects 
at  these  structures. 

Literature  Reviews  by  University  of  Houston 
Researchers 

The  objective  of  this  section  is  to  investigate  current 
debris-flow  management  strategies  that  can  be  applied  to  alle¬ 
viate  LWC  damage.  Flows  of  water  charged  with  substantial 
quantities  of  sediments  ranging  from  fine  clay  sizes  to  boulders 
are  termed  “debris  flows”  (Pierson  and  Costa,  1987;  Fiebiger, 
1997;  VanDine  and  others,  1997;  Wieczorek  and  others,  1997). 
Debris-flow  management  is  an  important  component  of  struc¬ 
tural  design  in  watersheds  with  substantial  bed  mobility.  Plan¬ 
ning  and  construction  of  countermeasures  have  evolved  from 
consolidation  and  retaining  functions  (large  debris  dams)  in  the 


early  1970s  to  more  elaborate  and  integrated  systems  that 
include  initiation  prevention,  magnitude  reduction,  and  flow 
diversion  (Fiebiger,  1997). 

Debris-Flow  Management  for  Low-Water  Crossing 
Protection 

A  debris-flow  management  system  can  be  divided  into 
prevention,  which  deals  with  watershed  management,  and  con¬ 
trol.  There  are  two  styles  of  control:  passive  and  active.  Passive 
control  includes  landscape  management,  forest  management, 
and  hazard  and  risk  mapping.  The  purpose  of  such  control  is  to 
avoid  the  effects  of  debris  flow  by  zoning  or  warning  (VanDine 
and  others,  1997).  Active  control,  which  tries  to  minimize  the 
areal  extent  of  the  debris  flow,  involves  stabilization  of  debris 
sources,  energy  dissipators,  dam  systems,  diversion  structures, 
or  discharge  control  features  (Fiebiger,  1997).  Active  measures 
also  are  classified  according  to  areas  of  application,  including 
the  source  or  initiation  area,  transportation  zone,  and  deposition 
zone  (Campbell,  1975;  Baldwin  and  others,  1987;  Hungr  and 
others,  1987).  Thus,  active  measures  within  source  areas 
include  landscape  and  forest  management  according  to  Hungr 
and  others  (1987).  Their  intended  function  is  reduction  of 
debris-flow  production.  Active  control  in  the  transportation 
zone  would  include  energy  dissipators  and  diversion  structures. 
Energy  dissipators  reduce  the  flow  magnitude  and  convey  the 
debris  surge  to  designated  area.  Finally,  dam  systems  and  the 
associated  discharge  and  runoff  control  structures,  which  are 
designed  to  control  the  areal  extent  of  the  debris  flow,  would  be 
classified  as  active  measures  applicable  in  the  deposition  zone. 
An  integrated  system  can  be  built  by  combining  several  of  the 
different  management  tools  listed  in  table  5  to  achieve  an  effec¬ 
tive  control  structure.  Selection  and  design  criteria  are  site 
specific,  taking  into  account  the  character  of  the  flow,  the  chan¬ 
nel,  and  the  debris  fan.  The  complexity  of  the  design  depends 
on  the  structure  or  entity  to  protect  and  on  the  available  capital 
(VanDine  and  others,  1997).  The  following  sections  discuss 
flow  management  methods  relevant  to  the  project  that  prompted 
this  investigation.  Watershed,  landscape,  and  forest  manage¬ 
ment  will  not  be  presented  because  their  implementation  is 
beyond  the  scope  of  the  project. 

Passive-control  hazard  mapping  and  zoning  starts  by  the 
recognition  of  destruction  patterns  of  debris  flows:  direct 
impact,  indirect  impact,  and  flood  zones  (Hungr  and  others, 
1987).  The  direct  impact  zone  is  delineated  by  high  energy  and 
destructive,  rapidly  moving,  high-discharge  pulses  that  influ¬ 
ence  the  deposition  area  just  downstream  of  transportation 
zones.  The  direct  impact  zone,  where  deposition  begins,  is 
located  in  areas  where  slopes  decrease  below  the  transportation 
critical  level,  which  is  very  site  specific.  Inspection  of  nearby 
deposits  that  are  evidence  of  debris  flow  should  give  an  indica¬ 
tion  of  the  critical  slope.  In  selecting  that  slope,  attention  should 
be  paid  to  channel  confinement,  which  can  be  more  critical  than 
the  slope  (Hungr  and  others,  1987).  Indirect  impacts,  including 
burial  of  structures,  result  from  deposition  of  large  debris 
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Table  5.  Debris-flow  management  system  diagram  (Hungr  and  others,  1987;  Fiebiger,  1997;  VanDine,  1997). 


Debris-flow  prevention 


Debris-flow  control 


Watershed  management 


Passive  control 


Active  control 


1.  Landscape  management  (source  zone) 

2.  Forest  management  (source  zone) 

•  reforestation/controlled  harvest 

•  forest  road  construction  control 

3.  Hazard  and  risk  mapping  (source,  transport,  and  deposition 
zone) 


1.  Debris  source  stabilization 

•  source  area  grading 

•  riprap 

•  baffles 

•  check  dams 

•  creek  channel  linings 

2.  Energy  dissipaters  (transportation  zone) 

•  debris-flow  breakers 

•  debris-flow  screens  or  rakes 


3.  Dam  systems  (deposition  zone) 

•  debris-flow  structure 

•  retention  dams 

•  steel  grid  dams 

•  open  debris  basin 

•  close  debris  basin 


4.  Diversion  structures  (transportation  zone) 

•  deflection  walls 

•  deflection  berms 

•  terminal  berms 


5.  Discharge  control  features 


volumes  moving  at  low  velocities  after  passing  the  direct 
impact  zone.  Those  deposits  can  develop  bars  and  temporary 
earth  dams  that  divert  the  flow  to  an  unexpected  area  and  create 
flood  zones.  Indirect  impacts  and  flood  zones  are  linked  to  the 
runout  length,  which  in  turn  is  related  to  the  volume  of  debris, 
magnitude  at  the  beginning  of  the  deposition  zone,  nature  of  the 
debris  flow,  and  slope  of  the  deposition  area.  Equations  to 
derive  deposition  thickness  and  runout  length  are  presented  by 
Innes  (1983),  Jeyapalan  and  others  (1983),  and  Hungr  and  oth¬ 
ers  (1984).  Once  the  hazard  map  is  created,  development  in  dan¬ 
gerous  areas  is  restricted  by  zoning  the  areas  as  hazardous,  park 
land,  or  agricultural  areas  (Hungr  and  others,  1987).  Mapping 
also  can  contribute  to  public  infrastructure  design  by  determin¬ 
ing  areas  of  high  sediment  loads. 

Active  control  consists  of  debris  source  stabilization, 
energy  dissipators,  dam  systems,  diversion  structures,  discharge 
control  features,  and  combinations  of  one  or  more  of  those 
structures. 

Debris  source  stabilization  is  important.  Steep  slopes  are 
one  of  the  control  elements  in  the  initiation  and  amplification  of 
debris-flow  magnitude.  Grading  the  source  area  to  a  more  uni¬ 
form  and  smooth  slope  is  a  common  practice  to  mitigate  sliding 
hazard.  However,  the  soil  cementing  vegetation  requires  seed¬ 
ing  in  spring  or  early  summer  to  allow  full  growth  before  the 
arrival  of  heavy  rains  to  prevent  severe  erosion  (Baldwin  and 
others,  1987).  Riprap  also  can  be  used  on  over-steepened  slopes 


(fig.  7).  As  shown  in  the  figure,  the  riprap  needs  to  be  keyed  into 
the  bedrock  at  a  recommended  depth  of  about  1  meter  to  avoid 
sliding.  When  the  finished  slope  is  greater  than  50  percent, 
Baldwin  and  others  (1987)  recommend  covering  the  surface 
with  chain-link  fence. 

Scour  of  loose  blocks  of  sediment  and  slides  of  channel 
banks  are  an  important  supply  of  debris  flow  (Hungr  and  others, 
1987;  Fiebiger,  1997).  Scour  and  slides  can  be  reduced  by  ele¬ 
vating  the  water  surface  level.  This  is  achieved  by  building 
check  dams  along  steep  channels.  The  method  is  very  efficient 
but  quite  expensive.  As  alternatives,  steel  cable  nets  (DeNatale, 
1997)  and  steel  bags  of  cobbles  (Baldwin  and  others,  1987) 
often  are  used. 

Baffles  can  be  built  in  the  transportation  zone  to  deflect, 
check,  or  regulate  debris  flows.  Baffles  usually  are  made  of 
treated  timber  or  steel.  Baffles  also  have  been  used  as  temporary 
energy  dissipators,  designed  to  fail  at  the  impact  while  offering 
some  resistance  to  the  incoming  flow.  In  this  case,  baffles  are 
set  up  in  a  sand-filled  foundation  rather  than  in  concrete  (Bald¬ 
win  and  others,  1987). 

Channel  lining  is  used  to  stabilize  the  bed,  protect  the 
banks,  and  maintain  channel  alignment.  These  stabilizations  are 
important  for  LWC  abutment  protection.  The  lining  can  be 
made  of  riprap  that  subsequently  is  grouted  with  concrete. 

To  limit  lining  damage  from  debris  slides  and  bank  slumps, 
excessive  velocities  generated  by  the  smoothness  of  concrete. 
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Figure  7.  Riprap  (left)  on  hillside  in  Pacifica,  California,  and  riprap  design  recommendation  (right)  (from  Baldwin  and  others,  1987). 


and  maintenance  because  of  abrasion,  masonry  linings  are 
used  in  Europe  and  boulders  embedded  in  steel  fiber-reinforced 
concrete  are  used  in  British  Columbia  (Hungr  and  others,  1987). 

Energy  dissipaters  decrease  the  energy  of  the  incoming 
flow  by  means  of  flow  breakers  that  can  stand  alone  (fig.  8)  or 
can  be  coupled  with  other  structures  to  enhance  efficiency  or 
accomplish  other  controls.  When  built  as  stand-alone  structures, 
flow  breakers  are  made  of  massive  concrete  piers  raised  at  a 
design  height  across  the  channel  bed.  The  combination  of  flow 
breakers  with  debris  dams  dramatically  reduces  debris-flow 
rates  and  thus  decreases  the  size  needed  for  the  structure.  The 
combination  also  can  serve  as  debris  traps  to  retain  oversize 
material.  Low-rise  flow  breakers  also  can  be  built  at  the  bottom 


of  settling  basins  to  serve  as  deposit  structures.  Sediment  can  be 
removed  after  each  flood. 

Dam  systems  shown  in  figure  9  can  have  a  complex  struc¬ 
ture  or  consist  of  a  simple  excavation  and  piled  up  berms  around 
the  selected  site,  known  as  debris  dams.  Debris  dams  often  are 
coupled  with  other  debris-flow  control  structures,  such  as  flow 
breakers  and  rakes  or  screens,  to  accomplish  several  functions, 
such  as  stabilizing  the  channel  bed  and  decreasing  slope  ero¬ 
sion,  by  elevating  the  water  surface  level,  sorting  oversize  mate¬ 
rial,  and  retaining  debris  overload. 

Open  debris  basins  are  designed  with  a  large  surface 
area  to  maximize  deposition  area  and  minimize  depth  so  that 
some  downstream  sediment  supply  still  occurs.  The  basins 


Figure  8.  Stand-alone  flow  breakers  (left)  and  flow  breakers  in  debris  dam  (right)  (from  Fiebiger,  1997). 
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Figure  9.  Debris-flow  dam  with  a  screen  of  steel  beams  (left,  from  Fiebiger,  1997)  and  open  debris  basin  design  (right,  from  VanDine, 
1997). 


are  built  across  the  existing  channel  path  when  diversion  is  not 
necessary. 

Closed  debris  basins  (fig.  10)  have  a  more  complex 


deflection  berms,  which  would  canalize  debris  flows  toward 
them.  To  protect  the  community  of  Lions  Bay,  British 
Columbia,  from  debris  flows  in  Harvey  Creek,  a  closed  debris 


structure.  In  addition  to  the  excavated  basin  and  the  containing  basin  was  built  in  the  mid-1980s.  The  basin  can  contain 


berms,  they  have  debris  straining  features  that  are  built  on-line. 
The  debris  straining  features  can  be  raised  at  the  end  run  of 


2,648,600  cubic  feet  of  debris  and  the  dam  is  49  feet  high.  It  is 
an  earth-filled  structure  with  a  complex  and  heavy  reinforced 
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Figure  10.  Closed  debris  basin  schematic  (left)  and  closed  debris  basin  on  Harvey  Creek  in  British  Columbia  (right)  (from  VanDine,  1997). 
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Figure  11.  Deflection  berm  schematic  (left)  and  deflection  berm  on  Boulder  Creek  in  British  Columbia  (right)  (from  VanDine,  1997). 


concrete  face  that  allows  passage  of  the  normal  creek  and 
strained  debris  flows. 

Diversion  structures  are  composed  of  deflection  and  termi¬ 
nal  berms  and  open  and  closed  debris  basins  (V anDine  and 
others,  1997).  Deflection  berms  (fig.  1 1)  are  built  to  deflect  the 
flow  path  while  providing  lateral  boundaries  to  avoid  areal 
spreading  of  debris.  The  figure  also  shows  a  deflection  berm 
built  on  Boulder  Creek,  British  Columbia  (Hawley,  1989; 
VanDine  and  others,  1997),  to  prevent  debris  flows  from 
overtopping  the  highway.  The  berm  was  10  to  56  feet  high 
and  1 ,640  feet  long,  the  lower  end  composed  of  riprap.  The 


structure  was  designed  to  handle  about  530,000  cubic  feet  of 
debris. 

Steel  wire  bags  of  cobbles  can  be  used  as  deflection 
devices  (fig.  12).  The  bags  are  effective  in  sustaining  high 
impact  forces.  The  bags  are  deformed  during  high-magnitude 
debris  flows.  However  these  deformations  do  not  affect  bag 
efficiency.  The  bags  are  easily  replaced  and  can  be  reused 
after  maintenance.  As  shown  in  figure  12,  high  segments  are 
designed  to  sustain  impacts  and  decrease  potential  for  over¬ 
topping.  The  lower  segments  should  be  embedded  into  the 
bedrock  to  avoid  sliding. 


-GRADED  PAD- 


Figure  12.  Gabion  deflection  wall,  Tiburon,  California  (left),  and  design  gabion  deflection  wall  (right)  (from  Baldwin  and  others,  1987). 
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Figure  13.  Terminal  berm  schematic  (left)  and  terminal  berm  on  Cypre  Creek  in  British  Columbia  (right)  (from  VanDine,  1997). 


Terminal  berms  (fig.  13)  are  earth-filled  debris  trap 
belts  built  as  far  down  as  possible  from  the  avalanche  scarp 
to  maximize  deposition  area  and  decrease  impact  force.  The 
berms  are  built  on-line  of  the  existing  channel,  enclosing  a  dep¬ 
osition  area  big  enough  to  contain  the  expected  maximum  vol¬ 
ume  of  debris.  Weirs  are  made  on  the  sides  to  allow  the  creation 
of  secondary  channels  that  would  move  relatively  low  debris- 
loaded  flows  downstream.  Figure  13  shows  a  terminal  berm 
constructed  on  the  west  coast  of  Vancouver  Island  in  1989  on  a 
seasonal  tributary  to  Cypre  Creek.  It  is  656  feet  long  and 
11.5  feet  high  and  withstood  a  debris-flow  volume  of  176,573 
cubic  feet. 

Research  on  Debris-Flow  Countermeasures 

Flexible  wire  rope  barriers  have  been  used  successfully 
in  controlling  debris  flows.  Wire  ropes  were  used  in  Pacifica, 
California,  as  fences  in  baffle  systems,  which  were  designed 
to  reduce  flow  velocity  (Baldwin  and  others,  1987).  They 
also  served  as  debris  straining  structures  in  closed  debris  basins 
(VanDine  and  others,  1997).  In  1994,  a  flexible  wire  rope 
prevented  about  2,100  cubic  feet  of  debris  flows  from  causing 
damage  to  California  State  Route  41.  Following  the  incident, 
the  USGS  conducted  field  experiments  to  study  the  response 
of  flexible  wire  rope  barriers  to  debris-flow  charges  (DeNatale 
and  others,  1997).  A  reinforced  concrete  flume  tilted  at  a  5- 
percent  slope  was  installed  near  Blue  Forest,  Oregon.  The  chan¬ 
nel  was  6.5  feet  wide,  312  feet  long,  and  3.9  feet  deep.  A  run¬ 
out  pad,  82  feet  long  was  mounted  at  the  flume  exit.  The  pad 
was  tilted  at  a  5-percent  slope  (fig.  14).  The  system  could 
deliver  44  cubic  feet  of  debris  flow  moving  at  32.8  feet  per  sec¬ 


ond.  In  order  to  create  a  soil  slurry  that  would  be  similar  to 
debris-flow  composition  and  properties  and  that  would  fit 
within  the  flume  operational  constraints,  the  experimental 
soil  was  artificially  made  of  42.1-percent  gravel,  73.9-percent 
sand,  and  2.0-percent  silt  and  clay.  The  typical  grain-size 
diameter  was  6.2  millimeters  and  the  effective  grain  size  was 
0.35  millimeter. 

Four  models  of  barriers  were  tested.  The  first  two  models 
were  wire  rope  netting  8  by  30  feet  with  chain  links  to  create 
nets  of  12  by  12  inches  and  8  by  8  inches  (fig.  15).  Ropes  of 
0.75-inch  diameter  were  used  for  the  frame  and  the  columns 
were  W4  by  13-inch  steel.  The  latter  two  models  were  both  8  by 
30  inches  with  a  0.75-inch  diameter  rope  frame  and  W8  by  48- 
inch  steel  support  columns.  The  6-by  6-inch  clear  opening  of 
the  third  model  was  made  with  chicken  wire.  The  fourth  model 
had  12-inch  opening  interlocking  rings  and  a  silt  screen  as  an 
additional  feature.  Flow  velocities  were  measured  using  1- 
square-meter  gridlines  marked  on  the  run-out  pad,  a  time- 
stamped  videotape  recorder,  and  an  overhead  camera.  Move¬ 
ment  of  markers,  including  sprayed  gravels  and  chunks,  were 
recorded.  Barrier  deformations  were  measured  with  a  cable 
extensometer  attached  to  the  nets.  After  running  multiple  debris 
flows  into  the  barriers,  only  models  3  and  4  did  not  suffer  col¬ 
umn  or  anchor-cable  connection  failure.  Deformations  were 
less  than  6.6  feet  and  only  0.05  percent  of  debris  passed  under¬ 
neath  the  barriers.  The  investigators  (DeNatale  and  others, 
1997)  recommended  that  wire  rope  barriers  could  be  used  to 
mitigate  small  debris  flows. 

Flow  breakers  are  energy  dissipaters  that  decrease  the 
energy  of  the  incoming  flow.  They  can  be  used  as  stand-alone 
structures  or  be  coupled  with  other  control  measures.  Chen  and 
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Figure  14.  U.S.  Geological  Surveyflume  on  Blue  River  at  Blue  Forest,  Oregon  (left),  and  schematic  of  the  experiment  to  study  the 
response  of  flexible  wire  rope  barriers  to  debris-flow  charges  (right)  (from  DeNatale  and  others,  1997). 


Ho  (1997)  studied  effects  of  four  different  types  of  flow 
breakers  that  the  authors  referred  to  as  “dams”  (gravitational, 
grid,  step-grid,  and  column)  on  debris  flows.  As  shown  in 
figure  16,  a  gravitational  flow  breaker  is  a  heavy  reinforced 
concrete  structure.  The  model  used  for  experiments  had  a 
5:3  front  slope,  clear  openings,  and  heights  1.5  and  5  times 
greater  than  the  largest  grain  of  the  experimental  debris  flow. 
The  largest  grain  diameter  was  38  millimeters.  Grid  and 
step-grid  flow  breakers  were  piles  of  steel-bar  cubes  with  open¬ 
ings  1.5  times  greater  than  the  largest  particle  size.  Different 
numbers  of  layers,  or  square  units  making  up  a  cube,  were 


tested.  The  column  flow  breaker  was  made  of  equally  spaced 
steel  bars  of  a  length  3  times  greater  than  the  maximum  diame¬ 
ter  size. 

A  slurry,  used  as  debris  flow,  was  prepared  and  compacted 
to  a  concentration  of  about  1 19  pounds  (mass)  per  foot  with 
4-percent  water  and  placed  at  the  top  of  a  flume  6.7  feet  long, 
tilted  at  34  percent  (fig.  17).  After  compaction,  the  slurry  was 
saturated  with  water  for  24  hours  (VanDine,  1986;  Baldwin  and 
others,  1987).  The  compacted  slurry  then  was  flushed  down¬ 
ward  as  water  rushed  in  the  flume  at  about  49  cubic  feet  per 
second. 


Figure  15.  Flexible  wire  rope  barriers  with  netting  (left;  models  1, 2,  and  3)  and  with  interlocking  rings  (right;  model  4)  (from  DeNatale 
and  others,  1997). 
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Figure  16.  Gravitational  flow  broaker  (top  loft),  five-layer  grid  flow  breaker  (top  right),  step-grid  flow  breaker  (bottom  left),  and  column 
flow  breaker  (bottom  right)  (from  Chen  and  Ho,  1997). 


Results  showed  that  the  more  layers  (cubes)  the  grid  flow 
breakers  were  made  of,  the  more  debris  could  be  retained.  How¬ 
ever  there  is  a  number  of  layers  above  which  retaining  effi¬ 
ciency  could  not  be  improved.  The  column  debris-flow  breaker 
was  the  most  efficient  in  decreasing  debris-flow  energy  but  did 
not  retain  debris.  Experiments  also  revealed  that  gravitational 
flow  breakers  were  best  suited  as  a  structure  outlet.  Fiebiger 
(1997)  had  the  same  conclusion  about  gravitational  flow 
breakers. 


Lin  and  others  (1997)  studied  retaining  rates  and  segrega¬ 
tion  efficiency  of  three  types  of  flow  breakers  (fig.  18) — the 
column,  which  Chen  and  Ho  (1997)  also  studied,  the  hori¬ 
zontal  column  referred  to  as  a  beam,  and  a  more  elaborated 
version  of  the  step-grid  Chen  and  Ho  (1997)  studied.  In 
Chen  and  Ho  (1997)  experiments,  column  and  grid  flow 
breakers  have  clear  openings  3  and  1.5  times  as  large  as  the 
largest  grain-size  diameter,  respectively.  In  Lin  and  others 
(1997),  the  openings  of  the  horizontal  column  and  grid 


Figure  17.  Flume  test  apparatus  (from  Chen  and  Ho,  1997). 
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Figure  18. 


Column  flow  broakor  (left),  beam  flow  breaker  (center),  and  grid  flow  breaker  (right)  (from  Lin  and  others,  1997). 


flow  breakers  are  2  times  greater  than  the  largest  grain-size 
diameter. 

A  debris-flow  slurry  was  made  by  compaction  of  the  sedi¬ 
ment,  submersion  for  2  hours,  then  drainage.  Successive  tests 
were  run  using  sediments  with  maximum  grain-size  diameters 
of  25.4  and  50.8  millimeters.  The  slurry  was  placed  in  the  upper 
part  of  a  flume  tilted  at  23,  27,  31,  and  then  35  percent.  Debris 
flows  were  triggered  with  water  flushed  in  at  flow  rates  of  0.07 
and  0.14  cubic  feet  per  second,  successively. 

The  retaining  rate  ( R )  was  defined  as  the  ratio  of  weight 
of  sand  retained  to  the  sum  of  sand  passed  and  sand  retained. 
Segregation  efficiency  was  quantified  by  the  separation  index, 
defined  as  the  ratio  of  average  grain  size  of  sand  that  passed  to 
average  grain  size  of  the  original  sample.  Experiments  showed 
through  regression  analysis  that  retaining  rate  was  highly  corre¬ 
lated  to  flume  tilt  as  was  the  ratio  of  flow  breaker  opening  to 
maximum  grain  size  ( L/Dmax );  whereas,  separation  index  was 
well  correlated  to  L/Dmax  .  Other  important  conclusions  are 
that  pore  pressure  is  controlled  by  seepage,  not  static  water 
pressure,  and  that  impact  force  magnitude  increases  drastically 
when  coarse  sediment  strikes  the  flow  breakers.  The  authors  of 
this  report  also  noted  that  recorded  dynamic  impact  forces  were 
more  than  2  times  as  much  as  theoretical  values. 

Mizuyama  and  Mizuno  (1997)  noticed  that  flow  breakers 
were  used  more  often  because  of  their  ability  to  allow  for  a  suf¬ 
ficient  sediment  supply  downstream  as  well  as  for  fish  to  pass 
through  the  openings.  However,  because  material  of  apprecia¬ 
ble  size  can  pass  through  flow  breakers,  prediction  of  the  flow 
hydrograph  and  sediment  content  is  necessary.  The  authors  set 
up  an  experiment  to  compare  measured  instantaneous  sediment 
discharge  to  a  theoretically  derived  passing  hydrograph.  The 
peak  sediment  reduction  rate  (P)  is  defined  as: 

P  =  \-R,  (20) 

where  R  is  the  retaining  rate  (see  Lin  and  others,  1997)  defined 
as: 
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where  Qsp  is  the  peak  sediment  flow  without  structure,  and 
Qspo  is  the  peak  sediment  flow  with  structure.  From  previous 
studies,  Mizuyama  and  Mizuno  (1997)  noted  that  the  reduction 
rate  (P )  was  a  function  of  L/Dmax  (see  Lin  and  others,  1997) 
and  sediment  concentration  ( C )  and  then  established  an  empir¬ 
ical  equation  relating  these  parameters, 
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where  P  =  1  -  Qsp/Qspo  ■  Assuming  that  the  flow  was  not 
influenced  by  deposits  inside  the  flume  and  that  trapped  sedi¬ 
ment  never  flowed  out,  Mizuyama  and  Mizuno  (1997)  estab¬ 
lished  the  following  equation: 
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where  Qs(t)  is  the  predicted  instantaneous  debris-flow 
discharge  through  a  structure,  and  Qso(t )  is  the  instantaneous 
debris-flow  discharge  without  a  structure.  Dmax  is  taken  as 
Dg5  .  For  experimental  verification,  the  authors  (this  report)  set 
up  a  wooden  flume  that  was  13  feet  long  and  6  inches  wide, 
tilted  at  17,  27,  and  36  percent.  Mizuyama  and  Mizuno  (1997) 
measured  Qso(t)  and  Qs-measured(t) ;  the  authors  (this  report) 
subsequently  computed  Qs-computeii(t)  and  find  good 
agreement  between  computed  and  measured  debris-flow 
hydrographs. 


Conclusions 

Numerous  debris-flow  mitigation  strategies  have  been 
proposed  around  the  world.  The  strategies  span  from  simple  pil¬ 
ing  of  berms  used  as  low-cost  flow  deflectors  or  debris  trap 
basins  to  more  complex  structures  that  exceed  millions  of 
dollars. 

Current  research  focuses  on  flow  breakers  that  can  be  used 
as  stand-alone  structures  or  as  discharge  control  structures. 
USGS  experiments  with  the  possibility  of  using  flexible  rope 
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wires  as  barriers  showed  that  6-  by  6-inch  opening  wire  rope 
netting  with  chicken  wire  and  interlocking-ring  wire  ropes  were 
promising  but  should  be  limited  to  mitigation  of  small  debris 
flows.  Column  debris-flow  breakers  are  used  more  often 
because  they  stand  out  as  the  most  efficient  in  decreasing 
debris-flow  energy  and  allow  both  a  sufficient  sediment  supply 
downstream  and  for  fish  to  pass  through  the  openings.  Experi¬ 
ments  also  showed  that  retaining  rates  of  column  debris-flow 
breakers  were  highly  related  to  hill  slope,  the  ratio  of  flow 
breaker  opening  to  maximum  grain  size  (L/Dmax ),  and 
sediment  concentration.  Material  of  considerable  size  can  pass 
through  flow  breakers  and  affect  downstream  structures.  To 
have  better  control  over  the  inflow  and  outflow  of  debris  move¬ 
ments,  an  attempt  to  predict  the  movement  of  debris  flows  past 
flow  breakers  was  made  and  the  results  were  encouraging. 

Debris-flow  magnitude  should  be  alleviated  with  mitiga¬ 
tion  structures  (discussed  herein)  before  reaching  LWCs.  Coun¬ 
termeasures  at  LWCs  should  be  limited  to  erosion  and  structural 
stability  controls.  Important  considerations  in  designing  debris- 
flow  mitigation  structures  would  be  to  quantify  impact  forces, 
the  deposit  volume  per  event,  and  the  magnitude  of  discharges, 
frequency  of  events,  and  sediment  size  distribution  of  debris 
flows. 

Literature  Reviews  by  U.S.  Geological  Survey 
Researchers 

Edwards  and  Glysson  (1999)  serves  as  the  guide  for 
USGS  investigations  involving  the  sampling  of  both  suspended 
and  bed  sediment.  Equipment  description  and  technique  discus¬ 
sions  are  included.  This  report  is  a  comprehensive  methodolog¬ 
ical  guide.  The  authors  distinguish  bed-material  sampling  from 
bed-load  sampling,  the  former  associated  with  bed  composition 
and  the  latter  defined  as  “sediment  that  moves  by  sliding,  roll¬ 
ing,  or  bouncing  along  on  or  near  the  streambed.’’ 

Edwards  and  Glysson  (1999)  report  that  current  equipment 
for  bed-material  sampling  can  accurately  represent  particles 
finer  than  16  millimeters  in  diameter  and  only  can  accommo¬ 
date  particles  finer  than  30  to  40  millimeters.  A  variety  of  equip¬ 
ment  is  available  for  this  purpose,  including  hand-held  samplers 
and  cable-and-reel  samplers.  Standard  equipment  for  sampling 
bed  material  comprising  gravel,  cobbles,  and  boulders  is  not 
available,  and  apparatus  designed  for  this  purpose  are  experi¬ 
mental.  For  particles  greater  than  16  millimeters,  the  popular 
and  randomized  Wolman  pebble  count  (Wolman,  1954)  is  rec¬ 
ommended.  The  pebble  count  is  conducted  by  establishing  a 
grid  of  sample  sites,  after  which  100  pebbles  are  chosen  ran¬ 
domly  at  each  site.  The  long,  intermediate,  and  short  axes  of 
each  clast  are  measured  and  recorded,  the  clast  then  is  placed 
into  a  size  interval,  and  the  percentage  of  each  interval  provides 
a  frequency  histogram  and  a  grain- size  curve.  Another  method 
mentioned  involves  the  use  of  a  device  known  as  a  “Zeiss  par¬ 
ticle-size  analyzer”  in  combination  with  photography  of  the 
streambed.  Finally,  the  used  of  a  “pipe  dredge”  is  mentioned  for 
non-wadable  streams. 


Edwards  and  Glysson  (1999)  report  that  bed-load 
sampling  is  difficult,  as  a  device  on  or  near  the  bed  disturbs 
velocity  distribution  and  rate  of  bed-load  movement.  Compli¬ 
cating  this  is  the  spatial  and  temporal  variation  of  both  the  flow 
velocity  near  the  bed  and  the  bed  load.  For  best  results,  a  bed¬ 
load  sampler  must  sample  the  mass  or  volume  of  particles 
through  a  given  width  in  a  specified  period  of  time.  Direct- 
collecting  samplers  fall  within  four  categories:  (1)  box  or  bas¬ 
ket,  (2)  pan  or  tray,  (3)  pressure  difference,  and  (4)  slot  or  pit 
samplers.  The  efficiency  of  bed-load-sampling  devices  is 
defined  as  the  “ratio  of  the  mass  of  bed  load  collected  to  the 
mass  of  the  bed  load  that  would  have  passed  through  the  sam¬ 
pler  width  in  the  same  time  period  had  the  sampler  not  been 
there.”  For  the  greatest  efficiency,  the  slot  or  pit  sampler  is  pre¬ 
ferred,  however,  extraction  of  the  sample  can  be  either  difficult 
or  costly.  Pressure-difference  samplers  are  hydrodynamically 
designed  to  maintain  ambient  entrance  flow  velocities  by  creat¬ 
ing  a  pressure  drop  at  the  exit.  The  Helley-Smith  bed-load 
sampler  is  the  most  popular  type  of  pressure-difference  sam¬ 
pler.  It  consists  of  an  expanding  3-  by  3-inch  nozzle,  sample 
bag,  and  frame,  allowing  for  the  collection  of  clasts  less  than  76 
millimeters  at  mean  velocities  less  than  9.8  feet  per  second.  An 
enlarged  6-  by  6-inch  version  is  available  for  streams  with 
larger  particle  sizes.  These  samplers  are  cumbersome  and 
require  a  cable-reel  suspension  system.  Helley-Smith  bed-load 
samplers  are  categorized  by  the  ratio  of  the  entrance  nozzle 
expansion,  the  original  sampler  having  a  3.22  area  ratio.  The 
USGS  endorses  the  use  of  a  Helley-Smith  bed-load  sampler 
with  a  1.40  area  ratio  of  the  expanding  entrance  nozzle, 
although  data  obtained  with  the  original  3.22  area  ratio  will  still 
be  accepted.  Efficiencies  have  been  determined  as  about  1.0, 
however,  it  should  be  noted  that  the  sampler  has  the  potential  to 
collect  some  particles  moving  in  suspension. 

Edwards  and  Glysson  (1999)  emphasize  the  temporal  vari¬ 
ability  of  bed-load  transport,  associated  with  its  characteristic 
movement  in  the  form  of  dunes,  cycles,  or  slugs.  Additionally, 
the  cross-channel  variation  is  typified  by  very  little  movement 
near  the  banks  to  maximum  transport  toward  mid-stream  or  in 
the  thalweg.  To  obtain  a  reasonable  estimate  of  the  mean  cross- 
sectional  bed-load  discharge,  the  collection  of  about  40  rate 
measurements  is  recommended,  subdividing  the  cross  section 
into  20  verticals  and  making  two  passes.  The  time  that  the  sam¬ 
pler  rests  on  the  bed  should  be  equal  for  all  verticals  in  a  cross 
section  but  does  not  have  to  be  the  same  for  both  passes.  The 
sampling  time  is  preferred  to  be  no  more  than  60  seconds,  but  it 
should  be  short  enough  to  fill  no  more  than  40  percent  of  the 
sample  bag  from  the  vertical  with  the  highest  transport  rate. 
This  is  referred  to  as  the  single  equal-width-increment  (SEWI) 
method.  The  multiple  equal-width-increment  (MEWI)  method 
collects  one  sample  at  four  or  more  evenly  spaced  verticals, 
followed  by  repeated  passes  until  a  total  of  40  samples  are 
collected.  Finally,  the  unequal-width-increment  (UWI)  method 
collects  one  sample  at  four  or  more  unevenly  spaced  verticals, 
followed  by  repeated  passes  until  a  total  of  40  samples  are 
collected. 
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Edwards  and  Glysson  (1999)  report  that  the  bed-load 
transport  rate  at  a  vertical  is: 

M, 

R,  =  K—  ,  (24) 

ti 

where  R ,  is  the  bed-load  transport  rate  at  vertical  i ,  in  tons  per 
day  per  foot;  M,  is  mass  of  sample  collected  at  vertical  i ,  in 
grams;  tl  is  the  time  the  sampler  was  on  the  bottom  at  vertical 
i ,  in  seconds;  and  K  is  a  unit  conversion  factor  used  to  convert 
grams  per  second  per  foot  into  tons  per  day  per  foot.  For  a  3- 
inch  nozzle  K  =  0.381 ,  and  for  a  6-inch  nozzle  K  =  0.190. 

The  cross-sectional  bed-load  discharge  is  calculated  using 
the  total  cross-section  method,  if  the  sample  times  at  each  ver¬ 
tical  were  equal,  the  verticals  were  evenly  spaced,  and  the  first 
sample  was  collected  at  one-half  the  sample  width  from  the 
bank.  If  any  of  these  conditions  are  not  met,  either  the  midsec¬ 
tion  or  the  mean-section  method  should  be  used.  The  total 
cross-section  formula  is: 

W 

Qb  =  K— Mt,  (25) 

‘t 

where  QB  is  the  bed-load  discharge,  in  tons  per  day;  WT  is  total 
width  of  stream,  in  feet;  tT  is  the  total  time  that  the  sampler  is 
on  the  bed,  in  seconds;  MT  is  the  total  mass  of  sample  collected 
from  all  verticals  samples  in  the  cross  section,  in  grams;  and  K 
is  defined  as  in  eq.  24. 

Stevens  and  Yang  (1989)  provide  a  summary  of  fluvial 
sediment  discharge  equations,  compare  measured  and  com¬ 
puted  results  from  the  equations,  and  provide  guidance  to  select 
equations  for  a  variety  of  flow  and  sediment  conditions.  Sedi¬ 
ment  discharge  equations  fall  into  three  categories:  determinis¬ 
tic,  probabilistic,  and  regression.  The  deterministic  approach 
assumes  a  correlation  between  dependent  and  independent  vari¬ 
ables  and  is  advantageous  if  the  independent  variables  are 
known.  Probability-based  equations  rely  on  predictions  of  par¬ 
ticle  motion,  and  regression  analysis  is  used  to  derive  empirical 
relations  between  sediment  discharge  rates  and  some  other 
parameter.  The  authors  include  13  fluvial  sediment-discharge 
equations,  five  of  which  are  bed-load  discharge  and  eight  are 
bed-material  discharge.  The  bed-load  discharge  is  the  sediment 
moving  in  continuous  contact  with  the  bed.  The  bed-material 
discharge  is  the  sediment  readily  exchangeable  with  the  bed- 
material,  moving  both  as  bed  load  and  in  suspension. 

Stevens  and  Yang  (1989),  through  their  review  of  sedi¬ 
ment  discharge  equations  accompanied  by  measured  values, 
computed  discrepancy  ratios  between  equations  and  made 
recommendations  regarding  use  of  the  equations.  A  discrep¬ 
ancy  ratio,  an  indicator  of  the  equation  accuracy,  is  the  ratio 
between  computed  and  measured  sediment  discharges.  Stevens 
and  Yang  (1989)  conclude  that  reliable  bed-load-discharge 
equations  are  available  in  Schoklitsch  (1934)  and  Meyer-Peter 
and  Muller  (1948).  Meyer-Peter  and  Muller  (1948)  is  recom¬ 
mended  for  bed  material  coarser  than  5  millimeters.  The 


Einstein  (1950)  equation  is  best  suited  for  bed-load-dominated 
rivers,  and  the  Yang  (1984)  equation  can  be  applied  when  most 
bed  materials  range  from  2  to  10  millimeters.  Stevens  and  Yang 
(1989)  conclude  that  reliable  bed-material-discharge  equations 
are  available  in  Yang  (1973)  (sand  only),  Engelund  and  Hansen 
(1967),  Ackers  and  White  (1973),  and  Yang  (1984)  (gravel 
only). 

Elliott  (2002)  addresses  a  problem  notably  similar  to 
that  of  TxDOT  project  0-4695.  In  this  study,  the  entrainment 
potential  is  computed  for  a  river  having  a  frequently  mobile 
streambed  composed  of  gravel,  cobbles,  and  boulders.  Cross- 
sectional  surveys  and  water-surface  elevations  were  used  to 
compute  the  boundary  shear  stress  at  different  locations,  includ¬ 
ing  the  mean  streambed,  thalweg,  and  bar  surfaces.  The  equa¬ 
tion  for  boundary  shear  stress  is: 

X0  =  JDS,  (26) 

where  xo  is  mean  boundary  shear  stress,  in  newtons  per  square 
meter;  y  is  specific  weight  of  water  (9,807  newtons  per  square 
meter);  D  is  the  mean  depth  of  flow,  in  meters;  and  S  is  the 
energy  gradient,  in  meters  per  meter.  The  slope  of  the  water  sur¬ 
face  often  is  used  as  a  substitute  for  S . 

Elliott  (2002)  then  uses  the  critical  shear  stress  to  describe 
the  point  at  which  entrainment  of  bed  material  is  initiated.  The 
critical  shear  stress  is  calculated  with  the  following  equation: 

U  =  Veils -y)d50,  (27) 

where  xc  is  the  critical  shear  stress,  in  newtons  per  square 
meter;  x*  is  the  dimensionless  critical  shear  stress  or  Shields 
parameter  (Shields,  1936;  Neill,  1968);  y  is  the  specific  weight 
of  sediment  (often  taken  as  2.65  times  the  specific  weight  of 
water);  y  is  the  specific  weight  of  water  (9,807  newtons  per 
square  meter);  and  d50  is  the  median  particle  size,  in  meters. 

Elliott  (2002)  concludes  that  sediment-entrainment  poten¬ 
tial  for  any  given  geomorphic  surface  or  flow  can  be  expressed 
as  the  ratio  of  boundary  shear  stress  to  critical  shear  stress 
( T0/xc ).  The  condition  at  which  the  median  particle  size 
becomes  entrained  is  referred  to  as  the  partial  entrainment 
threshold  (xo/xc  =  1 ).  The  complete  entrainment  threshold 
( xfl/xc  =  2 )  describes  the  point  at  which  complete  or  wide¬ 
spread  mobilization  of  the  median  particle  size  occurs. 

Gustavson  (1978)  is  the  only  identified  citation  on  gravel 
entrainment  in  the  Edwards  Plateau.  Using  historic  streamflow 
data  of  the  Nueces  River  in  Texas,  specifically  for  large  floods, 
the  author  computed  shear  stresses  to  estimate  stream  compe¬ 
tence.  The  results  show  that  the  magnitude  of  flow  required  to 
initiate  movement  of  gravel  with  an  intermediate  diameter  of 
about  0.5  inch  (about  1.2  centimeters)  occurs  only  0.1  percent 
of  the  time  each  year,  and  a  12-percent  probability  exists  each 
year  for  floods  capable  of  transporting  (initiating  movement  of) 
almost  all  of  the  bed  material.  Additionally,  the  river  can  trans¬ 
port  over  one-half  of  its  bed  material  during  flows  with  a  rela¬ 
tively  frequent  recurrence  interval  (1-  to  2-year  annual  peak 
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Figure  19.  Processed  La ndsat  Thematic  Mapping  scene  for  part  of  western  Arizona  (from  Mayer  and  Pearthree,  2002,  pi.  3).  The  image 
shows  changes  caused  by  Hurricane  Nora  in  different  colors.  The  colors  result  from  processing  bands  4  and  5  before  and  afterthe  storm. 
The  yellow  color  highlights  areas  in  which  flooding  deposited  fine-grained  material.  Mayer  and  Pearthree  (2002)  conclude  that  green  and 
cyan  colors  probably  highlight  areas  of  vegetation  growth. 


discharge)  and  can  mobilize  more  than  90  percent  of  its  bed 
material  once  every  7  or  8  years. 

Biedenharn  and  others  (2000)  provide  practical  guidance 
for  calculation  of  the  effective  discharge  of  rivers.  Effective  dis¬ 
charge  is  defined  as  the  streamflow  able  to  transport  the  most 
sediment,  specifically  the  bed-load  material.  Using  data  avail¬ 
able  from  gaging  stations,  a  standardized  procedure  was  devel¬ 
oped  to  calculate  effective  discharge.  In  general,  the  approach 
uses  a  flow-duration  curve  and  a  sediment-transport  rating 
curve.  For  ungaged  sites,  nearby  gaging  stations  within  the 
same  drainage  basin  can  be  used  or  a  regionalized  flow-duration 
curve  must  be  developed.  Following  the  creation  of  a  flow- 
duration  curve  from  measured  discharge,  a  sediment-transport 
equation  is  used  or  a  sediment-transport  rating  curve  must 
be  created  through  intensive  data  collection.  The  sediment- 
transport  rating  curve  has  sediment-load  data  (vertical  axis) 
plotted  against  discharge  (horizontal  axis).  A  power  function 
through  regression  often  provides  the  sediment-transport  equa¬ 
tion.  Finally,  a  bed-material  load  histogram  is  created  to  repre¬ 
sent  the  total  bed-material  load  transported  by  each  discharge 
class.  This  histogram  is  produced  by  multiplying  the  bed- 
material  load  by  the  frequency  of  occurrence  of  its  correspond¬ 
ing  discharge.  The  effective  discharge  corresponds  to  the  modal 
class  of  the  histogram. 


Mayer  and  Pearthree  (2002)  describe  the  use  of  standard 
Landsat  Thematic  Mapping  (TM)  data  to  detect  sediment 
changes  on  the  land  surface  caused  by  flooding  in  the  deserts  of 
the  southwestern  United  States.  Mayer  and  Pearthree  (2002) 
also  provide  additional  references  for  space-based  flood  map¬ 
ping  and  resultant  changes  to  landscape.  The  authors  process 
Landsat  TM  scenes  recorded  in  July  1997  and  in  October  1997 
over  the  Tiger  Wash  area  of  southwestern  Arizona  (fig.  19). 
This  area  experienced  widespread  flooding  as  a  result  of 
Hurricane  Nora  in  late  September  1997.  Landsat  TM  data  is 
available  from  the  USGS  National  Center  for  Earth  Resources 
Observational  Science  (EROS).  The  Landsat  satellites  are  about 
700  kilometers  above  the  earth  and  revisit  the  same  area  every 
16  days. 

Mayer  and  Pearthree  (2002)  process  infrared  bands  4  and 
5.  Each  band  has  a  spatial  resolution  of  30  meters.  The  authors 
processed  the  scenes  to  accentuate  scene  changes  that  are 
believed  to  reflect  changes  to  sediment  and  vegetation.  The  pro¬ 
cessing  is  complex  and  includes  adjustments  for  such  factors  as 
inherent  differences  in  solar  illumination  and  geographic  refer¬ 
encing  errors.  Extensive  field  analysis  was  used  to  map  regions 
of  sediment  change  following  the  Hurricane  Nora  flooding 
event.  Although  the  meaning  of  the  changes  in  Landsat  TM 
scenes  on  a  pixel  to  pixel  basis  are  clear,  the  changes  between 
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scenes  are  “strikingly”  similar  to  the  field  mapping  of  sediment 
changes.  Mayer  and  Pearthree  (2002,  pi.  3)  provide  an  spectac¬ 
ular  high-resolution  color  plate  of  a  part  of  western  Arizona 
showing  the  sediment  and  other  changes  caused  by  Hurricane 
Nora  flooding.  They  conclude  that  Landsat  TM  data  with  spe¬ 
cific  processing  can  be  used  to  detect  regional  changes  to  chan¬ 
nels  and  overbank  areas  resulting  from  substantial  flooding  in 
ephemeral  streams  in  the  southwestern  United  States. 

Mayer  and  Pearthree  (2002)  is  a  promising  reference  that 
describes  application  of  geographic  information  systems  (GIS) 
and  remotely  sensed  data  suitable  to  study  the  bed-load- 
transport  phenomena  in  the  Edwards  Plateau  area  in  Texas. 
Because  of  the  substantial  rearrangement  of  the  bed  during 
large  floods  in  the  area,  new  insights  in  the  phenomena  could 
be  made  through  periodic  processing  of  Landsat  data  or  high- 
resolution  digital  orthophoto  quarter-quadrangles  (DOQQs). 
However,  a  specific  methods  and  processing  study  of  Landsat 
TM  data  over  the  Edwards  Plateau  in  the  context  of  bed-load 
transport  is  required  before  cost-effective  periodic  data  process¬ 
ing  is  feasible,  DOQQs  might  be  the  most  cost-effective,  high- 
resolution  imagery  to  use  in  spatial  and  temporal  analyses  of 
bed-load-transport. 

Summary 

A  review  of  the  literature  addressing  sediment  transport  in 
gravel-bed  river  systems  and  structures  designed  to  control  bed¬ 
load  mobility  is  provided  as  part  of  Texas  Department  of  Trans¬ 
portation  research  project  0-4695:  Guidance  for  Design  in 
Areas  of  Extreme  Bed-Load  Mobility.  The  study  area  com¬ 
prises  Edwards,  Kimble,  and  Real  Counties  in  the  Edwards 
Plateau  in  central  Texas.  The  primary  focus  of  the  literature 
review  is  on  journal  articles,  edited  volumes,  and  government 
publications.  The  literature  review  provides  an  outline  and 
foundation  for  the  research  project  to  characterize  extreme  bed¬ 
load  mobility  in  rivers  and  streams  across  the  study  area.  The 
literature  review  also  provides  a  basis  for  potential  modifica¬ 
tions  to  low-water  stream-crossing  design  in  the  study  area. 

The  literature  associated  with  sediment  transport  in  gravel- 
bed  rivers  and  application  of  theory  to  river  management  is 
extensive,  global,  and  multi-disciplinary.  Investigations  focus 
on  data  derived  from  natural  channels  or  experimental  flumes, 
and  approaches  vary  from  deterministic  to  empirical.  A  variety 
of  methods  are  used  to  study  gravel-bed  river  characteristics, 
including  surveying  channel  geometry,  particle  tracing,  bed- 
material  and  bed-load  sampling,  and  developing  sediment- 
discharge  rating  curves,  among  others.  Additionally,  efforts  to 
manipulate  the  transport  of  sediment  in  gravel-bed  rivers 
include  the  construction  of  debris  traps,  velocity  control  struc¬ 
tures,  or  flow  diversion  structures.  Regardless  of  the  method 
implemented  to  alter  the  sediment  regime  of  gravel-bed  rivers, 
a  thorough  knowledge  of  watershed  characteristics,  local  pro¬ 
cesses  governing  the  transport,  and  collection  of  data  character¬ 


izing  sediment  mobility  is  essential  for  the  success  of  any  man¬ 
agement  practice. 

An  assessment  of  the  literature  provides  a  guide  for 
research  of  sediment-transport  problems  in  gravel-bed  rivers. 
Problems  associated  with  extreme  bed-load  mobility  in  rivers 
and  streams  of  the  Edwards  Plateau  would  benefit  from  a  robust 
geographic,  field,  and  laboratory  assessment  of  the  phenome¬ 
non.  Geographic  techniques  include  geographic  information 
systems  (GIS)  analysis  of  field  data,  aerial  photographic  inter¬ 
pretation  and  time-series  analysis  of  sediment  condition  at 
selected  sites,  and  measurement  of  channel  and  basin  character¬ 
istics.  Field  procedures  should  include  bed-material  particle- 
size  analyses;  surveys  of  channel  geometry,  bed  configuration, 
and  slope;  bed-load  measurements  accompanied  by  flow  veloc¬ 
ity  values;  and  particle  tracing.  Especially  important  are  mea¬ 
surements  during  high-flow  events  and  post-flood  comparisons 
to  earlier  conditions.  Field  measurements  could  be  used  to 
develop  sediment-transport  rating  curves,  calculate  shear 
stresses  and  stream  power,  compare  with  channel  geometry  or 
bed  reconfiguration,  and  facilitate  the  design  of  structures  or 
improvements  to  stream  crossings.  In  comparison  with  field 
assessments,  laboratory  flume  investigations  have  the  advan¬ 
tage  of  controlling  desired  parameters,  such  as  particle-size  or 
channel  width,  to  test  the  response  of  other  variables.  Addition¬ 
ally,  flume  experiments  involving  structural  models  might  be 
conducted  to  test  the  effectiveness  of  possible  designs  or 
improvements  to  stream  crossings. 
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•  The  authors  describe  the  spatial  and  temporal  definition 

of  stability. 
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permits  sedimentation. 

Kjartanson,  B.H.,  Heilers,  G.A.,  Lohnes,  R.A.,  and  Klaiber, 
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